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Abstract 

In the first part of this contribution, we review the development of the theory of 
scale relativity and its geometric framework constructed in terms of a fractal and 
nondifferentiable continuous space-time. This theory leads (i) to a generalization 
of possible physically relevant fractal laws, written as partial differential equation 
acting in the space of scales, and (ii) to a new geometric foundation of quantum 
mechanics and gauge field theories and their possible generalisations. 

In the second part, we discuss some examples of application of the theory to 
various sciences, in particular in cases when the theoretical predictions have been 
validated by new or updated observational and experimental data. This includes 
predictions in physics and cosmology (value of the QCD coupling and of the cosmo- 
logical constant), to astrophysics and gravitational structure formation (distances of 
extrasolar planets to their stars, of Kuiper belt objects, value of solar and solar-like 
star cycles), to sciences of life (log-periodic law for species punctuated evolution, 
human development and society evolution) , to Earth sciences (log-periodic decelera- 
tion of the rate of California earthquakes and of Sichuan earthquake replicas, critical 
law for the arctic sea ice extent) and tentative applications to system biology. 

1 Introduction 



One of the main concern of the theory of scale relativity is about the foundation of 
quantum mechanics. As it is now well known, the principle of relativity (of motion) 
underlies the foundation of most of classical physics. Now, quantum mechanics, though it 
is harmoniously combined with special relativity in the framework of relativistic quantum 
mechanics and quantum field theories, seems, up to now, to be founded on different 
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grounds. Actually, its present foundation is mainly axiomatic, i.e., it is based on postulates 
and rules which are not derived from any underlying more fundamental principle. 

The theory of scale relativity suggests an original solution to this fundamental problem. 
Namely, in its framework quantum mechanics may indeed be founded on the principle of 
relativity itself, provided this principle (applied up to now to position, orientation and 
motion) be extended to scales. One generalizes the definition of reference systems by 
including variables characterizing their scale, then one generalizes the possible transfor- 
mations of these reference systems by adding, to the relative transformations already 
accounted for (translation, velocity and acceleration of the origin, rotation of the axes), 
the transformations of these variables, namely, their relative dilations and contractions. 
In the framework of such a newly generalized relativity theory, the laws of physics may 
be given a general form that transcends and includes both the classical and the quantum 
laws, allowing in particular to study in a renewed way the poorly understood nature of 
the classical to quantum transition. 

A related important concern of the theory is the question of the geometry of space-time 
at all scales. In analogy with Einstein's construction of general relativity of motion, which 
is based on the generalization of flat space-times to curved Riemannian geometry, it is 
suggested, in the framework of scale relativity, that a new generalization of the description 
of space-time is now needed, toward a still continuous but now nondifferentiable and frac- 
tal geometry (i.e., explicitly dependent on the scale of observation or measurement). New 
mathematical and physical tools are therefore developed in order to implement such a gen- 
eralized description, which goes far beyond the standard view of differentiable manifolds. 
One writes the equations of motion in such a space-time as geodesies equations, under the 
constraint of the principle of relativity of all scales in nature. To this purpose, covariant 
derivatives are constructed that implement the various effects of the nondifferentiable and 
fractal geometry. 

As a first theoretical step, the laws of scale transformation that describe the new 
dependence on resolutions of physical quantities are obtained as solutions of differential 
equations acting in the space of scales. This leads to several possible levels of description 
for these laws, from the simplest scale invariant laws to generalized laws with variable 
fractal dimensions, including log-periodic laws and log-Lorentz laws of "special scale- 
relativity", in which the Planck scale is identified with a minimal, unreachable scale, 
invariant under scale transformations (in analogy with the special relativity of motion in 
which the velocity c is invariant under motion transformations). 

The second theoretical step amounts to describe the effects induced by the internal 
fractal structures of geodesies on motion in standard space (of positions and instants). 
Their main consequence is the transformation of classical dynamics into a generalized, 
quantum-like self-organized dynamics. The theory allows one to define and derive from 
relativistic first principles both the mathematical and physical quantum tools (complex, 
spinor, bispinor, then multiplet wave functions) and the equations of which these wave 
functions are solutions: a Schrodinger-type equation (more generally a Pauli equation 
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for spinors) is derived as an integral of the geodesic equation in a fractal space, then 
Klein-Gordon and Dirac equations in the case of a full fractal space-time. We then briefly 
recall that gauge fields and gauge charges can also be constructed from a geometric re- 
interpretation of gauge transformations as scale transformations in fractal space-time. 

In a second part of this review, we consider some applications of the theory to various 
sciences, particularly relevant to the questions of evolution and development. In the 
realm of physics and cosmology, we compare the various theoretical predictions obtained 
at the beginning of the 90's for the QCD coupling constant and for the cosmological 
constant to their present experimental and observational measurements. In astrophysics, 
we discuss applications to the formation of gravitational structures over many scales, with 
a special emphasis on the formation of planetary systems and on the validations, on the 
new extrasolar planetary systems and on Solar System Kuiper belt bodies discovered since 
15 years, of the theoretical predictions of scale relativity (made before their discovery). 
This is completed by a vahdation of the theoretical prediction obtained some years ago 
for the solar cycle of 11 yrs on other solar-like stars whose cycles are now measured. In 
the realm of life sciences, we discuss possible applications of this extended framework to 
the processes of morphogenesis and the emergence of prokaryotic and eukaryotic cellular 
structures, then to the study of species evolution, society evolution, embryogenesis and 
cell confinement. This is completed by applications in Earth sciences, in particular to 
a prediction of the Arctic ice rate of melting and to possible predictivity in earthquake 
statistical studies. 

2 Theory 

2.1 Foundations of scale relativity theory 

The theory of scale relativity is based on the giving up of the hypothesis of manifold 
differentiability. In this framework, the coordinate transformations are continuous but 
can be nondifferentiable. This implies several consequences [69], leading to the following 
steps of construction of the theory: 

(1) One can prove the following theorem [69, 72, 7, 22, 23]: a continuous and nondif- 
ferentiable curve is fractal in a general meaning, namely, its length is explicitly dependent 
on a scale variable e, i.e., £ = C{e), and it diverges, £ — > oo, when e — 0. This theorem 
can be readily extended to a continuous and nondifferentiable manifold, which is therefore 
fractal, not as an hypothesis, but as a consequence of the giving up of an hypothesis (that 
of differentiability). 

(2) The fractality of space-time [69, 105, 66, 67] involves the scale dependence of 
the reference frames. One therefore adds to the usual variables defining the coordinate 
system, new variables e characterizing its 'state of scale'. In particular, the coordinates 
themselves become functions of these scale variables, i.e., X = X{e). 
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(3) The scale variables e can never be defined in an absolute way, but only in a relative 
way. Namely, only their ratio p = e' /e does have a physical meaning. In experimental 
situations, these scales variables amount to the resolution of the measurement apparatus 
(it may be defined as standard errors, intervals, pixel size, etc.). In a theoretical analysis, 
they are the space and time differential elements themselves. This universal behavior leads 
to extend the principle of relativity in such a way that it applies also to the transformations 
(dilations and contractions) of these resolution variables [67, 68, 69]. 



2.2 Laws of scale transformation 

2.2.1 Fractal coordinate and differential dilation operator 

Consider a variable length measured on a fractal curve, and, more generally, a non- 
differentiable (fractal) curvilinear coordinate £(s, e), that depends on some parameter s 
which characterizes the position on the curve (it may be, e.g., a time coordinate), and on 
the resolution e. Such a coordinate generalizes to nondifferentiable and fractal space-times 
the concept of curvilinear coordinates introduced for curved Riemannian space-times in 
Einstein's general relativity [69]. 

Such a scale-dependent fractal length £(s, e), remains finite and differentiable when 
£ 7^ 0, namely, one can define a slope for any resolution £, being aware that this slope is 
itself a scale-dependent fractal function. It is only at the limit e ^ that the length is 
infinite and the slope undefined, i.e., that nondifferentiability manifests itself. 

Therefore the laws of dependence of this length upon position and scale may be written 
in terms of a double differential calculus, i.e., it can be the solution of differential equations 
involving the derivatives of C with respect to both s and e. 

As a preliminary step, one needs to establish the relevant form of the scale variables 
and the way they intervene in scale differential equations. For this purpose, let us apply 
an infinitesimal dilation dp to the resolution, which is therefore transformed as e — e' = 
e{l + dp). The dependence on position is omitted at this stage in order to simplify the 
notation. By applying this transformation to a fractal coordinate £, one obtains, to first 
order in the differential element, 

C{e') = C{e + e dp) = C{e) + e dp = {1 + D dp) C{e), (1) 

where D is, by definition, the dilation operator. 

Since de/e = dine, the identification of the two last members of equation (1) yields 

This form of the infinitesimal dilation operator shows that the natural variable for the 
resolution is \ne, and that the expected new differential equations will indeed involve 



4 



quantities such as dC{s,e) / dine. This theoretical result agrees and explains the current 
knowledge according to which most measurement devices (of light, sound, etc.), including 
their physiological counterparts (eye, ear, etc..) respond according to the logarithm of 
the intensity (e.g., magnitudes, decibels, etc.). 



2.2.2 Self-similar fractals as solutions of a first order scale differential equa- 
tion 

Let us start by writing the simplest possible differential equation of scale, then by solving 
it. We shall subsequently verify that the solutions obtained comply with the principle 
of relativity. As we shall see, this very simple approach already yields a fundamental 
result: it gives a foundation and an understanding from first principles for self-similar 
fractal laws, which have been shown by Mandelbrot and many others to be a general 
description of a large number of natural phenomena, in particular biological ones (see, 
e.g., [60, 104, 59], other volumes of these series and references therein). In addition, the 
obtained laws, which combine fractal and scale-independent behaviours, are the equivalent 
for scales of what inertial laws are for motion [60]. Since they serve as a fundamental 
basis of description for all the subsequent theoretical constructions, we shall now describe 
their derivation in detail. 

The simplest differential equation of explicit scale dependence which one can write is 
of first order and states that the variation of C under an infinitesimal scale transformation 
dine depends only on C itself. Basing ourselves on the previous derivation of the form of 
the dilation operator, we thus write 

'-^-m. (3) 

The function /3 is a priori unknown. However, still looking for the simplest form of 
such an equation, we expand [3{C) in powers of £, namely we write /?(£) = a + hC + .... 
Disregarding for the moment the s dependence, we obtain, to the first order, the following 
linear equation, in which a and h are constants: 

a + hC. (4) 



dine 



In order to find the solution of this equation, let us change the names of the constants as 
tf = —b and Cq = a/rp, so that a + hC = —tf{C — Co). We obtain the equation 



dC 



Its solution reads 



-Tp dine. (5) 



m = Co {l+{-T\, (6) 
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where A is an integration constant. This solution corresponds to a length measured on a 
fractal curve up to a given point. One can now generalize it to a variable length that also 
depends on the position characterized by the parameter s. One obtains 



C{s,e) = Cois) 1+C(s) 




(7) 



in which, in the most general case, the exponent Tp may itself be a variable depending on 
the position. 

The same kind of result is obtained for the projections on a given axis of such a fractal 
length [69]. Let X{s,e) be one of these projections, it reads 



In this case Cx{s) becomes a highly fluctuating function which may be described by a 
stochastic variable. 

The important point here and for what follows is that the solution obtained is the 
sum of two terms, a classical-like, "differentiable part" and a nondifferentiable "fractal 
part", which is explicitly scale-dependent and tends to infinity when e [69, 17]. 
By differentiating these two parts in the above projection, we obtain the differential 
formulation of this essential result. 



where dx is a classical differential element, while is a differential element of fractional 
order. This relation plays a fundamental role in the subsequent developments of the 
theory. 

Consider the case when Tp is constant. In the asymptotic small scale regime, e <C A, 
one obtains a power-law dependence on resolution that reads 



We recognize in this expression the standard form of a self-similar fractal behaviour with 
constant fractal dimension Dp = 1 + Tp, which have already been found to yield a fair 
description of many physical and biological systems [60] . Here the topological dimension 
is Dp = 1, since we deal with a length, but this can be easily generalized to surfaces 
[Dp = 2), volumes [Dp = 3), etc., according to the general relation Dp = Dp + Tp. The 
new feature here is that this result has been derived from a theoretical analysis based on 
first principles, instead of being postulated or deduced from a fit of observational data. 

It should be noted that in the above expressions, the resolution is a length interval, 
e = SX defined along the fractal curve (or one of its projected coordinate). But one may 





dX = dx + d^, 



(9) 




(10) 
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also travel on the curve and measure its length on constant time intervals, then change 
the time scale. In this case the resolution e is a time interval, e = 6t. Since they are 
related by the fundamental relation 



SX^^ ~ 6t, (11) 
the fractal length depends on the time resolution as 



T 



1-1/Df 



X{s,6t)=Xo{s)xi^-j . (12) 

An example of the use of such a relation is Feynman's result according to which the mean 
square value of the velocity of a quantum mechanical particle is proportional to 6t~^ [33, 
p. 176], which corresponds to a fractal dimension Dp = 2, as later recovered by Abbott 
and Wise [1] by using a space resolution. 

More generally, (in the usual case when e = 6X), following Mandelbrot, the scale 
exponent Tp = Dp — Dp can be defined as the slope of the {lne,\nC) curve, namely 

ciln£ 

For a self-similar fractal such as that described by the fractal part of the above solution, 
this definition yields a constant value which is the exponent in Eq. (10). However, one can 
anticipate on the following, and use this definition to compute an "effective" or "local" 
fractal dimension, now variable, from the complete solution that includes the differentiable 
and the nondifferentiable parts, and therefore a transition to effective scale independence. 
Differentiating the logarithm of Eq. (16) yields an effective exponent given by 

The effective fractal dimension Dp = 1 + Tp therefore jumps from the nonfractal value 
Dp = Dp = 1 to its constant asymptotic value at the transition scale A. 



2.2.3 Galilean relativity of scales 

We can now check that the fractal part of such a law is compatible with the principle 
of relativity extended to scale transformations of the resolutions (i.e., with the principle 
of scale relativity). It reads C = Co{X/ey^ (Eq. 10), and it is therefore a law involving 
two variables {\nC and Tp) in function of one parameter (e) which, according to the 
relativistic view, characterizes the state of scale of the system (its relativity is apparent 
in the fact that we need another scale A to define it by their ratio). More generally, all 
the following statements remain true for the complete scale law including the transition 
to scale- independence, by making the replacement of £ by £ — Cq. Note that, to be 
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complete, we anticipate on what follows and consider a priori Tp to be a variable, even 
if, in the simple law first considered here, it takes a constant value. 

Let us take the logarithm of Eq. (10). It yields ln(£/£o) = tf ln(A/£:). The two 
quantities ln£ and Tp then transform, under a finite scale transformation e ^ e' = pe, 
as 

In— — = ln— Ti^lnp, (15) 

''-0 ''-0 

and, to be complete, 

r'p = Tp. (16) 

These transformations have exactly the same mathematical structure as the Galilean 
group of motion transformation (applied here to scale rather than motion), which reads 

x' = x-tv, t' = t. (17) 

This is confirmed by the dilation composition law, e —* e' ^ e", which writes 

e" e' e" 
In — = In - + In — , (18) 
e e e' 

and is therefore similar to the law of composition of velocities between three reference 
systems K, K' and K'" , 

V"{K"IK) = V{K'/K) + V'{K"/K'). (19) 

Since the Galileo group of motion transformations is known to be the simplest group that 
implements the principle of relativity, the same is true for scale transformations. 

It is important to realize that this is more than a simple analogy: the same physical 
problem is set in both cases, and is therefore solved under similar mathematical structures 
(since the logarithm transforms what would have been a multiplicative group into an 
additive group). Indeed, in both cases, it amounts to find the law of transformation of 
a position variable {X for motion in a Cartesian system of coordinates, hiC for scales 
in a fractal system of coordinates) under a change of the state of the coordinate system 
(change of velocity V for motion and of resolution In p for scale), knowing that these state 
variables are defined only in a relative way. Namely, V is the relative velocity between the 
reference systems K and K\ and p is the relative scale: note that e and e' have indeed 
disappeared in the transformation law, only their ratio remains. This remark founds the 
status of resolutions as (relative) "scale velocities" and of the scale exponent Tp as a "scale 
time" . 

Recall finally that, since the Galilean group of motion is only a limiting case of the 
more general Lorentz group, a similar generalization is expected in the case of scale 
transformations, which we shall briefly consider in Sec. 2.2.6. 
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2.2.4 Breaking of scale invariance 



The standard self-similar fractal laws can be derived from the scale relativity approach. 
However, it is important to note that Eq. (16) provides us with another fundamental 
result. Namely, it also contains a spontaneous breaking of the scale symmetry. Indeed, it 
is characterized by the existence of a transition from a fractal to a non-fractal behaviour 
at scales larger than some transition scale A. The existence of such a breaking of scale 
invariance is also a fundamental feature of many natural systems, which remains, in most 
cases, misunderstood. 

The advantage of the way it is derived here is that it appears as a natural, sponta- 
neous, but only effective symmetry breaking, since it does not affect the underlying scale 
symmetry. Indeed, the obtained solution is the sum of two terms, the scale-independent 
contribution ( different iable part), and the explicitly scale-dependent and divergent contri- 
bution (fractal part). At large scales the scaling part becomes dominated by the classical 
part, but it is still underlying even though it is hidden. There is therefore an apparent 
symmetry breaking, though the underlying scale symmetry actually remains unbroken. 

The origin of this transition is, once again, to be found in relativity (namely, in the 
relativity of position and motion). Indeed, if one starts from a strictly scale-invariant law 
without any transition, L = CQ{X/ey^ , then adds a translation in standard position space 
{C C + £i), one obtains 



Therefore one recovers the broken solution (that corresponds to the constant a 7^ in the 
initial scale differential equation). This solution is now asymptotically scale-dependent 
(in a scale-invariant way) only at small scales, and becomes independent of scale at large 
scales, beyond some relative transition Ai which is partly determined by the translation 
itself. 

2.2.5 Generalized scale laws 

Discrete scale invariance, complex dimension and log-periodic behaviour 

Fluctuations with respect to pure scale invariance are potentially important, namely the 
log-periodic correction to power laws that is provided, e.g., by complex exponents or 
complex fractal dimensions. It has been shown that such a behaviour provides a very 
satisfactory and possibly predictive model of the time evolution of many critical systems, 
including earthquakes and market crashes ([122] and references therein). More recently, 
it has been applied to the analysis of major event chronology of the evolutionary tree of 
life [19, 85, 86], of human development [14] and of the main economic crisis of western 
and precolumbian civilizations [44, 85, 50, 45]. 

One can recover log-periodic corrections to self-similar power laws through the re- 
quirement of covariance (i.e., of form invariance of equations) applied to scale differential 




(20) 
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equations [75]. Consider a scale-dependent function C{e). In the applications to temporal 
evolution quoted above, the scale variable is identified with the time interval |t — td, where 
tc is the date of a crisis. Assume that C satisfies a first order differential equation, 

-uC = 0, (21) 



dine 



whose solution is a pure power law C{e) oc e'^ (cf Sect. 2.2.2). Now looking for corrections 
to this law, one remarks that simply incorporating a complex value of the exponent u 
would lead to large log-periodic fiuctuations rather than to a controllable correction to 
the power law. So let us assume that the right-hand side of Eq. (21) actually differs from 
zero 

dC „ , , 

71 ^^ = X- 22 

dme 

We can now apply the scale covariance principle and require that the new function x 
be solution of an equation which keeps the same form as the initial equation 

u'x = 0. (23) 



dine 

Setting u' = u + 7] , we find that C must be solution of a second-order equation 

The solution reads C{e) = ae'^{l + be^), and finally, the choice of an imaginary exponent 
7] = iui yields a solution whose real part includes a log-periodic correction: 

C{e) = ae''[l + h cos(tu In e)] . (25) 

As previously recalled in Sect. 2.2.4, adding a constant term (a translation) provides a 
transition to scale independence at large scales. 

Lagrangian approach to scale laws In order to obtain physically relevant general- 
izations of the above simplest (scale- invariant) laws, a Lagrangian approach can be used 
in scale space, in analogy with its use to derive the laws of motion, leading to reverse the 
definition and meaning of the variables [75]. 

This reversal is an analog to that achieved by Galileo concerning motion laws. Indeed, 
from the Aristotle viewpoint, "time is the measure of motion". In the same way, the 
fractal dimension, in its standard (Mandelbrot's) acception, is defined from the topological 
measure of the fractal object (length of a curve, area of a surface, etc..) and resolution, 
namely (see Eq. 13) 

X ^ ^ d\\iC , . 

t=- ^ rp = DF-DT= , 26 

V dln{X/e) 
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In the case, mainly considered here, when L represents a length (i.e., more generally, a 
fractal coordinate), the topological dimension is = 1 so that Tp = Dp — 1. With 
Galileo, time becomes a primary variable, and the velocity is deduced from space and 
time, which are therefore treated on the same footing, in terms of a space-time (even 
though the Galilean space-time remains degenerate because of the implicitly assumed 
infinite velocity of light). 

In analogy, the scale exponent Tp = Dp — 1 becomes in this new representation a 
primary variable that plays, for scale laws, the same role as played by time in motion laws 
(it is called "djinn" in some publications which therefore introduce a five- dimensional 
'space-time-djinn' combining the four fractal fluctuations and the scale time). 

Carrying on the analogy, in the same way as the velocity is the derivative of position 
with respect to time, v = dx/dt, we expect the derivative of ln£ with respect to scale 
time Tp to be a "scale velocity". Consider as reference the self-similar case, that reads 
InC = Tp\n{\/e). Derivating with respect to Tp, now considered as a variable, yields 
dlnC/dTp = ln{X/e), i.e., the logarithm of resolution. By extension, one assumes that 
this scale velocity provides a new general definition of resolution even in more general 
situations, namely, 

' A \ dlnC 



One can now introduce a scale Lagrange function L(ln£, V, Tp), from which a scale action 
is constructed 



S= L{\nC,Y,Tp) dTp. (28) 

Jti 

The application of the action principle yields a scale Euler-Lagrange equation that writes 

dTpdY d\nC' ^ ' 

One can now verify that, in the free case, i.e., in the absence of any "scale force" (i.e., 
dL/dlnC = 0), one recovers the standard fractal laws derived hereabove. Indeed, in this 
case the Euler-Lagrange equation becomes 

dL/dY = const ^ V = const. (30) 

which is the equivalent for scale of what inertia is for motion. Still in analogy with motion 
laws, the simplest possible form for the Lagrange function is a quadratic dependence on 
the scale velocity, (i.e., L oc V^). The constancy of V = \n{X/e) means that it is 
independent of the scale time Tp. Equation (27) can therefore be integrated to give the 
usual power law behaviour, C = Co{\/eY^ , as expected. 

But this reversed viewpoint has also several advantages which allow a full implemen- 
tation of the principle of scale relativity: 
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(i) The scale time Tp is given the status of a fifth dimension and the logarithm of the 
resolution, V = \n{X/e), its status of scale velocity (see Eq. 27). This is in accordance with 
its scale-relativistic definition, in which it characterizes the state of scale of the reference 
system, in the same way as the velocity v = dx/dt characterizes its state of motion. 

(ii) This allows one to generalize the formalism to the case of four independent space- 
time resolutions, = \ii{X^/e^) = d\nC^ / drp. 

(iii) Scale laws more general than the simplest self-similar ones can be derived from 
more general scale Lagrangians [74, 75] involving "scale accelerations" W = d^lxiC/ drp = 
dln[X/e)/dTF, as we shall see in what follows. 

Note however that there is also a shortcoming in this approach. Contrarily to the 
case of motion laws, in which time is always fiowing toward the future (except possibly 
in elementary particle physics at very small time scales), the variation of the scale time 
may be non-monotonic, as exemplified by the previous case of log-periodicity. Therefore 
this Lagrangian approach is restricted to monotonous variations of the fractal dimension, 
or, more generally, to scale intervals on which it varies in a monotonous way. 



Scale dynamics The previous discussion indicates that the scale invariant behaviour 
corresponds to freedom (i.e. scale force-free behaviour) in the framework of a scale physics. 
However, in the same way as there are forces in nature that imply departure from iner- 
tial, rectilinear uniform motion, we expect most natural fractal systems to also present 
distorsions in their scale behaviour with respect to pure scale invariance. This implies 
taking non-linearity in the scale space into account. Such distorsions may be, as a first 
step, attributed to the effect of a dynamics of scale ("scale dynamics"), i.e., of a "scale 
field" , but it must be clear from the very beginning of the description that they are of ge- 
ometric nature (in analogy with the Newtonian interpretation of gravitation as the result 
of a force, which has later been understood from Einstein's general relativity theory as a 
manifestation of the curved geometry of space-time). 

In this case the Lagrange scale-equation takes the form of Newton's equation of dy- 
namics, 

dHnC , , 

^ = Pi) 

where /x is a "scale mass" , which measures how the system resists to the scale force, and 
where F = ci^ In £/(ir|. = ci ln(A/e) /dri? is the scale acceleration. 

In this framework one can therefore attempt to define generic, scale-dynamical be- 
haviours which could be common to very different systems, as corresponding to a given 
form of the scale force. 



Constant scale force A typical example is the case of a constant scale force. Setting 
G = F/fi, the potential reads cp = GlnC, in analogy with the potential of a constant 
force / in space, which is 9? = —fx, since the force is —dip/dx = f. The scale differential 
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equation writes 

It can be easily integrated. A first integration yields dlnC/drp = Grp + Vq, where Vq is 
a constant. Then a second integration yields a parabolic solution (which is the equivalent 
for scale laws of parabolic motion in a constant field), 

Y = Yo + Gtf ■ ln£ = ln/:o+Vor^ + ^Gr|, (33) 

where V = dlnC/dTp = \n{X/e). 

However the physical meaning of this result is not clear under this form. This is 
due to the fact that, while in the case of motion laws we search for the evolution of the 
system with time, in the case of scale laws we search for the dependence of the system 
on resolution, which is the directly measured observable. Since the reference scale A is 
arbitrary, the variables can be re-defined in such a way that Yq = 0, i.e., A = Aq. Indeed, 
from Eq. (33) one gets tf = {Y - Yo)/G = MX/e) - ln(A/Ao)]/G = \n{Xo/e)/G. Then 
one obtains 

1 / Ao \ , / £ \ 1 i_ 2 / ^0 



The scale time tf becomes a linear function of resolution (the same being true, as 
a consequence, of the fractal dimension Df = 1 + tf), and the (ln£,ln£:) relation is 
now parabolic instead of linear. Note that, as in previous cases, we have considered here 
only the small scale asymptotic behaviour, and that we can once again easily generalize 
this result by including a transition to scale-independence at large scale. This is simply 
achieved by replacing £ by (£ — Cq) in every equations. 

There are several physical situations where, after careful examination of the data, the 
power-law models were clearly rejected since no constant slope could be defined in the 
(log£, loge) plane. In the several cases where a clear curvature appears in this plane, e.g., 
turbulence [29], sandpiles [11], fractured surfaces in solid mechanics [13], the physics could 
come under such a scale-dynamical description. In these cases it might be of interest to 
identify and study the scale force responsible for the scale distorsion (i.e., for the deviation 
from standard scaling). 



2.2.6 Special scale-relativity 

Let us close this section about the derivation of scale laws of increasing complexity by 
coming back to the question of finding the general laws of scale transformations that meet 
the principle of scale relativity [68]. It has been shown in Sec. 2.2.3 that the standard 
self-similar fractal laws come under a Galilean group of scale transformations. However, 
the Galilean relativity group is known, for motion laws, to be only a degenerate form of 
the Lorentz group. It has been proved that a similar result holds for scale laws [68, 69]. 
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The problem of finding the laws of linear transformation of fields in a scale transfor- 
mation V = Inp (e — > e') amounts to finding four quantities, a(V), b{Y), c(V), and d{W), 
such that 

ln^ = a(V)ln^ + &(V)r^, (35) 

t'p = c(V) In ^ + d(y) Tp. 

Set in this way, it immediately appears that the current 'scale-invariant' scale trans- 
formation law of the standard form of constant fractal dimension (Eq. 15), given by 
a = l,6 = V, c = and d = 1, corresponds to a Galilean group. 

This is also clear from the law of composition of dilatations, e — >^ e' — > e", which has 
a simple additive form, 

V" = V + v. (36) 

However the general solution to the 'special relativity problem' (namely, find a, b, c and d 
from the principle of relativity) is the Lorentz group [58, 68]. This result has led to the 
suggestion of replacing the standard law of dilatation, e—*e' = gxehja. new Lorentzian 
relation, namely, for e < Xq and e' < Xq 

In — = ln(£/Ao) +lng 

""Ao l + lnf?ln(£/Ao)/ln2(A/Ao)' ^ ^ 

This relation introduces a fundamental length scale A, which is naturally identified, to- 
ward the small scales, with the Planck length (currently 1.6160(11) x 10~^^ m) [68], 

A = /p = {hG/cY^ (38) 

and toward the large scales (for e > Xq and e' > Xq) with the scale of the cosmological 
constant, L = A-^/^ [gg^ Qj^^p. 7.1]. 

As one can see from Eq. (37), if one starts from the scale e = A and apply any 
dilatation or contraction g, one obtains again the scale e' = A, whatever the initial value 
of Aq. In other words, A can be interpreted as a limiting lower (or upper) length-scale, 
impassable, invariant under dilatations and contractions. 

As concerns the length measured along a fractal coordinate which was previously scale- 
dependent as ln(£/£o) = lii(Ao/£:) for e < Aq, it becomes in the new framework, in the 
simplified case when one starts from the reference scale Cq 

In^ - , ^oMV^) (39) 



Co 



l-ln^(Ao/£)/ln^(Ao/A) 



The main new feature of scale relativity respectively to the previous fractal or scale- 
invariant approaches is that the scale exponent tf and the fractal dimension Dp = l + Tp, 
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which were previously constant {Dp = 2,tf = 1 ), are now exphcitly varying with scale, 
following the law (given once again in the simplified case when we start from the reference 
scale Co): 



Under this form, the scale covariance is explicit, since one keeps a power law form for the 
length variation, C = Co{\/ey''^^\ but now in terms of a variable fractal dimension. 

For a more complete development of special relativity, including its implications as 
regards new conservative quantities and applications in elementary particle physics and 
cosmology, see [68, 69, 72, 103]. 

The question of the nature of space-time geometry at Planck scale is a subject of intense 
work (see e.g. [5, 57] and references therein). This is a central question for practically all 
theoretical attempts, including noncommutative geometry [20, 21], supersymmetry and 
superstrings theories [43, 112], for which the compactification scale is close to the Planck 
scale, and particularly for the theory of quantum gravity. Indeed, the development of loop 
quantum gravity by Rovelli and Smolin [114] led to the conclusion that the Planck scale 
could be a quantized minimal scale in Nature, involving also a quantization of surfaces 
and volumes [115]. 

Over the last years, there has also been significant research effort aimed at the devel- 
opment of a 'Doubly-Special- Relativity' [4] (see a review in [5]), according to which the 
laws of physics involve a fundamental velocity scale c and a fundamental minimum length 
scale Lp, identified with the Planck length. 

The concept of a new relativity in which the Planck length-scale would become a 
minimum invariant length is exactly the founding idea of the special scale relativity theory 
[68], which has been incorporated in other attempts of extended relativity theories [15, 16]. 
But, despite the similarity of aim and analysis, the main difference between the 'Doubly- 
Special-Relativity' approach and the scale relativity one is that the question of defining 
an invariant length-scale is considered in the scale relativity/fractal space-time theory 
as coming under a relativity of scales. Therefore the new group to be constructed is a 
multiplicative group, that becomes additive only when working with the logarithms of 
scale ratios, which are definitely the physically relevant scale variables, as one can show 
by applying the Gell-Mann-Levy method to the construction of the dilation operator (see 
Sec. 2.2.1). 

2.3 Fractal space and quantum mechanics 

The first step in the construction of a theory of the quantum space-time from fractal 
and nondifferentiable geometry, which has been described in the previous sections, has 
consisted of finding the laws of explicit scale dependence at a given "point" or "instant" 
(under their new fractal definition). 




(40) 
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The next step, which will now be considered, amount to write the equation of motion in 
such a fractal space(-time) in terms of a geodesic equation. As we shall see, this equation 
takes, after integration the form of a Schrodinger equation (and of the Klein-Gordon and 
Dirac equations in the relativistic case). This result, first obtained in Ref. [69], has later 
been confirmed by many subsequent physical [72, 74, 28, 17] and mathematical works, in 
particular by Cresson and Ben Adda [22, 24, 8, 9] and Jumarie [51, 52, 53, 54], including 
attempts of generalizations using the tool of the fractional integro-differential calculus 
[9, 26, 54]. 

In what follows, we consider only the simplest case of fractal laws, namely, those char- 
acterized by a constant fractal dimension. The various generalized scale laws considered in 
the previous section lead to new possible generalizations of quantum mechanics [72, 103]. 

2.3.1 Critical fractal dimension 2 

Moreover, we simplify again the description by considering only the case Dp = 2. Indeed, 
the nondifferentiability and fractality of space implies that the paths are random walks 
of the Markovian type, which corresponds to such a fractal dimension. This choice is 
also justified by Feynman's result [33], according to which the typical paths of quantum 
particles (those which contribute mainly to the path integral) are nondifferentiable and of 
fractal dimension Dp = 2 [1]. The case Dp ^2, which yields generalizations to standard 
quantum mechanics has also been studied in detail (see [72, 103] and references therein). 
This study shows that Dp = 2 plays a critical role in the theory, since it suppresses the 
explicit scale dependence in the motion (Schrodinger) equation - but this dependence 
remains hidden and reappears through, e.g., the Heisenberg relations and the explicit 
dependence of measurement results on the reolution of the measurement apparatus. 

Let us start from the result of the previous section, according to which the solution 
of a first order scale differential equation reads for Dp = 2, after differentiation and 
reintro duct ion of the indices, 

dX" = dx" + rfe^ = v^'ds + C^/Kd's, (41) 

where Ac is a length scale which must be introduced for dimensional reasons and which, 
as we shall see, generalizes the Compton length. The C,^ are dimensionless highly fluctu- 
ating functions. Due to their highly erratic character, we can replace them by stochastic 
variables such that <C'^>= 0, <(C°)^>= -1 and <(C^)^>= 1 {k =1 to 3). The mean is 
taken here on a purely mathematic probability law which can be fully general, since the 
final result odes not depend on its choice. 

2.3.2 Metric of a fractal space-time 

Now one can also write the fractal fluctuations in terms of the coordinate differentials, 
d^^ = (^^vA'Nix^. The identification of this expression with that of Eq. (41) leads to 
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recover the Einstein-de Broglie length and time scales, 

X - _ _^ ^ _ K _ ^ 

^ dx/ds ' dt/ds E 

Let us now assume that the large scale (classical) behavior is given by Riemannian 
metric potentials g^^{x,y, z,t). The invariant proper time dS along a geodesic writes, in 
terms of the complete differential elements dX^^ = dx^ + d^^, 

dS^ = g^^dX^'dX'' = g^^idx^" + d^^'){dx'' + d^). (43) 

Now replacing the d^'s by their expression, one obtains a fractal metric [69, 84]. Its 
two-dimensional and diagonal expression, neglecting the terms of zero mean (in order to 
simplify its writing) reads 

dS^ = ^?oo(x, t) (l + Co' ^) c^de - g,,{x, t) [l + Cl ^ dx\ (44) 

We therefore obtain generalized fractal metric potentials which are divergent and ex- 
plicitly dependent on the coordinate differential elements [67, 69]. Another equivalent 
way to understand this metric consists in remarking that it is no longer only quadratic in 
the space-time differental elements, but that it also contains them in a linear way. 

As a consequence, the curvature is also explicitly scale-dependent and divergent when 
the scale intervals tend to zero. This property ensures the fundamentally non-Riemannian 
character of a fractal space-time, as well as the possibility to characterize it in an intrin- 
sic way. Indeed, such a characterization, which is a necessary condition for defining a 
space in a genuine way, can be easily made by measuring the curvature at smaller and 
smaller scales. While the curvature vanishes by definition toward the small scales in 
Gauss-Riemann geometry, a fractal space can be characterized from the interior by the 
verification of the divergence toward small scales of curvature, and therefore of physical 
quantities like energy and momentum. 

Now the expression of this divergence is nothing but the Heisenberg relations them- 
selves, which therefore acquire in this framework the status of a fundamental geometric 
test of the fractality of space-time [66, 67, 69]. 

2.3.3 Geodesies of a fractal space-time 

The next step in such a geometric approach consists in the identification of wave-particles 
with fractal space-time geodesies. Any measurement is interpreted as a selection of the 
geodesies bundle linked to the interaction with the measurement apparatus (that depends 
on its resolution) and/or to the information known about it (for example, the which- way- 
information in a two-slit experiment [72]. 

The three main consequences of nondifferentiability are: 
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(42) 



(i) The number of fractal geodesies is infinite. This leads to adopt a generalized 
statistical fluid-like description where the velocity V'^{s) is replaced by a scale-dependent 
velocity field V"^[X^(s, ds), s, ds]. 

(ii) There is a breaking of the reflexion invariance of the differential element ds. Indeed, 
in terms of fractal functions f{s,ds), two derivatives are defined, 

^ Xis + ds,ds)-X{s,ds)^ ^, ^ X{s,ds)-X{s-ds,ds)^ 

ds ds 

which transform one in the other under the reflection [ds ^ —ds), and which have a 
priori no reason to be equal. This leads to a fundamental two-valuedness of the velocity 
field. 

(iii) The geodesies are themselves fractal curves of fractal dimension Dp = 2 [33]. 
This means that one defines two divergent fractal velocity fields, V^[x{s,ds), s,ds] = 

v+[x{s),s] + w+[x{s,ds), s,ds] and V-[x{s,ds), s,ds] = V-[x{s),s] + W-[x{s,ds), s,ds], 
which can be decomposed in terms of differentiable parts v+ and V-, and of fractal parts 
w+ and W-. Note that, contrarily to other attempts such as Nelson's stochastic quantum 
mechanics which introduces forward and backward velocities [63] (and which has been 
later disproved [42, 126]), the two velocities are here both forward, since they do not 
correspond to a reversal of the time coordinate, but of the time differential element now 
considered as an independent variable. 

More generally, we define two differentiable parts of derivatives d+/ds and d-/ds, 
which, when they are applied to x^, yield the differential parts of the velocity fields, 
v1 = d^x^^/ds and = d^x^/ds. 

2.3.4 Covariant total derivative 

Let us first consider the non-relativistic case. It corresponds to a three-dimensional fractal 
space, without fractal time, in which the invariant ds is therefore identified with the time 
differential element dt. One describes the elementary displacements dX^, k = 1,2,3, on 
the geodesies of a nondifferentiable fractal space in terms of the sum of two terms (omitting 
the indices for simplicity) dX± = d±x + d^±, where dx represents the differentiable part 
and d^ the fractal (nondifferentiable) part, defined as 

d±x = v± dt, d^± = C± dt^^^. (46) 

Here (± are stochastic dimensionless variables such that <C±>= and <C±>= 1, and 
V is a. parameter that generalizes, up to the fundamental constant c/2, the Compton 
scale (namely, V = h/2m in the case of standard quantum mechanics). The two time 
derivatives are then combined in terms of a complex total time derivative operator [69], 



d 1 f d+ d^\ i / d^ d^ 
'dt~ 2\dt^ Ht) ~2\'dt~li 



(47) 



18 



Applying this operator to the differentiable part of the position vector yields a complex 
velocity 

In order to find the expression for the complex time derivative operator, let us first 
calculate the derivative of a scalar function /. Since the fractal dimension is 2, one needs 
to go to second order of expansion. For one variable it reads 

dt dt dX dt 2 • y ! 

The generalization of this writing to three dimensions is straighforward. 

Let us now take the stochastic mean of this expression, i.e., we take the mean on the 
stochastic variables C,± which appear in the definition of the fractal fiuctuation d^±. By 
definition, since dX = dx+d^ and <d^>= 0, we have <dX>= dx, so that the second term 
is reduced (in 3 dimensions) to v.Vf. Now concerning the term dX'^/dt, it is infinitesimal 
and therefore it would not be taken into account in the standard differentiable case. But in 
the nondifferentiable case considered here, the mean square fluctuation is non-vanishing 
and of order dt, namely, <d^'^>= 2T)dt, so that the last term of Eq. (49) amounts in 
three dimensions to a Laplacian operator. One obtains, respectively for the (+) and (-) 
processes, 

"^^^ ''^ +Vi.V±VA]f . (50) 



dt \dt 

Finally, by combining these two derivatives in terms of the complex derivative of Eq. (47), 
it reads [69] 

A = |. + v.V - zPA. (51) 
dt dt ^ ^ 

Under this form, this expression is not fully covariant [111], since it involves derivatives 

of the second order, so that its Leibniz rule is a linear combination of the flrst and second 

order Leibniz rules. By introducing the velocity operator [89] 

V = V - WV, (52) 

it may be given a fully covariant expression, 

namely, under this form it satisfles the flrst order Leibniz rule for partial derivatives. 

We shall now see that d/dt plays the role of a "covariant derivative operator" (in 
analogy with the covariant derivative of general relativity), namely, one may write in its 
terms the equation of physics in a nondifferentiable space under a strongly covariant form 
identical to the differentiable case. 
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2.3.5 Complex action and momentum 

The steps of construction of classical mechanics can now be followed, but in terms of 
complex and scale dependent quantities. One defines a Lagrange function that keeps its 
usual form, C{x,V,t), but which is now complex, then a generalized complex action 

S= r C{xy,t)dt. (54) 

Generalized Euler-Lagrange equations that keep their standard form in terms of the new 
complex variables can be derived from this action [69, 17], namely 

d dL dL 

From the homogeneity of space and Noether's theorem, one defines a generalized complex 
momentum given by the same form as in classical mechanics, namely. 

If the action is now considered as a function of the upper limit of integration in Eq. (54), 
the variation of the action from a trajectory to another nearby trajectory yields a gener- 
alization of another well-known relation of classical mechanics, 



V = V5. (57) 



2.3.6 Motion equation 



Consider, as an example, the case of a single particle in an external scalar field of poten- 
tial energy (but the method can be applied to any situation described by a Lagrange 
function). The Lagrange function , L = ^mv'^ — cj), is generalized as C{x, V, t) = ^mV^ — cj). 
The Euler-Lagrange equations then keep the form of Newton's fundamental equation of 
dynamics F = mdv/dt, namely, 

m-^V = -V0, (58) 

which is now written in terms of complex variables and complex operators. 

In the case when there is no external field (0 = 0), the covariance is explicit, since 
Eq. (58) takes the free form of the equation of inertial motion, i.e., of a geodesic equation, 

|V = 0, (59) 
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This is analog to Einstein's general relativity, where the equivalence principle leads to 
write the covariant equation of motion of a free particle under the form of an inertial 
motion (geodesic) equation Du^/ds = 0, in terms of the general-relativistic covariant 
derivative D, of the four- vector and of the proper time differential ds. 

The covariance induced by the effects of the nondifferentiable geometry leads to an 
analogous transformation of the equation of motions, which, as we show below, become af- 
ter integration the Schrodinger equation, which can therefore be considered as the integral 
of a geodesic equation in a fractal space. 

In the one-particle case the complex momentumP reads 

V = mV, (60) 

so that, from Eq. (57), the complex velocity V appears as a gradient, namely the gradient 
of the complex action 

V = VS/m. (61) 



2.3.7 Wave function 

Up to now the various concepts and variables used were of a classical type (space, 
geodesies, velocity fields), even if they were generalized to the fractal and nondifferen- 
tiable, explicitly scale-dependent case whose essence is fundamentally not classical. 

We shall now make essential changes of variable, that transform this apparently 
classical-like tool to quantum mechanical tools (without any hidden parameter or new 
degree of freedom). The complex wave function ip is introduced as simply another expres- 
sion for the complex action S, by making the transformation 

^|J = e'^/^\ (62) 

Note that, despite its apparent form, this expression involves a phase and a modulus since 
S is complex. The factor 5*0 has the dimension of an action (i.e., an angular momentum) 
and must be introduced because S is dimensioned while the phase should be dimensionless. 
When this formalism is applied to standard quantum mechanics, Sq is nothing but the 
fundamental constant h. As a consequence, since 

S = —iSo Inip, (63) 

one finds that the function ip is related to the complex velocity appearing in Eq. (61) as 
follows 

Q 

V = -i — VlnV^. (64) 
m 

This expression is the fondamental relation that connects the two description tools while 
giving the meaning of the wave function in the new framework. Namely, it is defined here 
as a velocity potential for the velocity field of the infinite family of geodesies of the fractal 
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space. Because of nondifferentiability, the set of geodesies that defines a 'particle' in this 
framework is fundamentally non-local. It can easily be generalized to a multiple particle 
situation, in particular to entangled states, which are described by a single wave function 
ip, from which the various velocity fields of the subsets of the geodesic bundle are derived 
as Vfc = —i (So/mk) Vk In?/;, where k is an index for each particle. The indistinguishability 
of identical particles naturally follows from the fact that the 'particles' are identified with 
the geodesies themselves, i.e., with an infinite ensemble of purely geometric curves. In 
this description there is no longer any point-mass with 'internal' properties which would 
follow a 'trajectory', since the various properties of the particle - energy, momentum, 
mass, spin, charge (see next sections) - can be derived from the geometric properties of 
the geodesic fluid itself. 



2.3.8 Correspondence principle 

Since we have V = —iSoVlmp = —iSo{'Vip)/ip, we obtain the equality [69] 

r^P = -ihVi) (65) 

in the standard quantum mechanical case Sq = h, which establishes a correspondence 
between the classical momentum p, which is the real part of the complex momentum in 
the classical limit, and the operator —iKV. 

This result is generalizable to other variables, in particular to the Hamiltonian. Indeed, 
a strongly covariant form of the Hamiltonian can be obtained by using the fully covariant 
form Eq. (53) of the covariant derivative operator. With this tool, the expression of the 
relation between the complex action and the complex Lagrange function reads 

Since V = VS and Ti = —dS/dt, one obtains for the generalized complex Hamilton 
function the same form it has in classical mechanics, namely [103, 96], 

n = V.V-C. (67) 

After expansion of the velocity operator, one obtains 7i = V.V — iW.V — C, which 
includes an additional term [111], whose origin is now understood as an expression of 
nondifferentiability and strong covariance. 



2.3.9 Schrodinger equation and Compton relation 

The next step of the construction amounts to write the fundamental equation of dynamics 
Eq. (58) in terms of the function ip. It takes the form 

2^o^(Vln^)=V0. (68) 
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As we shall now see, this equation can be integrated in a general way under the form of 
a Schrodinger equation. Replacing d/dt and V by their expressions yields 



V$ = iSq 



\ — (Vln^/^.V)(VlnV') +PA(VlnV') 
ot m 



This equation may be simplified thanks to the identity [69], 



(69) 



, , 2(Vln7/;.V)(VlnV^) + A(Vln7/^). (70) 
W J 

We recognize, in the right-hand side of Eq. (70), the two terms of Eq. (69), which were 
respectively in factor of S'o/m and V. This leads to definitely define the wave function as 

which means that the arbitrary parameter 5*0 (which is identified with the constant h in 
standard QM) is now linked to the fractal fluctuation parameter by the relation 

So = 2mV. (72) 

This relation (which can actually be proved instead of simply being set as a simplifying 
choice, see [94, 96]) is actually a generalization of the Compton relation, since the geo- 
metric parameter V =<dC,^> /2dt can be written in terms of a length scale as "D = Ac/2, 
so that, when Sq = h, it becomes A = h/mc. But a geometric meaning is now given 
to the Compton length (and therefore to the inertial mass of the particle) in the fractal 
space-time framework. 

The fundamental equation of dynamics now reads 



V0 = 2imV 



d 

— V \n^-i {2P(V In ^.V) (V In ^p) + PA(V In ^)} 



(73) 



Using the above remarkable identity and the fact that d/dt and V commute, it becomes 

-^ = -2VvU^\ni^ + V^\. (74) 
m y ot ijj ) 

The full equation becomes a gradient, 

v(l-2DvfiM£i±^U=o. (75) 



m \ J 

and it can be easily integrated, to finally obtain a generalized Schrodinger equation [69] 

V^A^ + tV^Jj-^^ = 0, (76) 
ot 2m 

up to an arbitrary phase factor which may be set to zero by a suitable choice of the 
phase. One recovers the standard Schrodinger equation of quantum mechanics for the 
particular case when V = h/2m. 
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2.3.10 Von Neumann's and Bern's postulates 



In the framework described here, "particles" are identified with the various geometric 
properties of fractal space (-time) geodesies. In such an interpretation, a measurement (and 
more generally any knowledge about the system) amounts to a selection of the sub-set of 
the geodesies family in which are kept only the geodesies having the geometric properties 
corresponding to the measurement result. Therefore, just after the measurement, the 
system is in the state given by the measurement result, which is precisely the von Neumann 
postulate of quantum mechanics. 

The Born postulate can also be inferred from the scale- relativity construction [17, 
94, 96]. Indeed, the probability for the particle to be found at a given position must 
be proportional to the density of the geodesies fiuid at this point. The velocity and the 
density of the fiuid are expected to be solutions of a Euler and continuity system of four 
equations, for four unknowns, {p,Vx,Vy,Vz). 

Now, by separating the real and imaginary parts of the Schrodinger equation, setting 
ijj = ^/P X e*^ and using a mixed representation {P,V), where V = {Vx,Vy,Vz}, one 
obtains precisely such a standard system of fiuid dynamics equations, namely. 



| + V-vV=-V - 21)^^ ' % + div(PV) = 0. (77) 



This allows one to univoquely identify P = with the probability density of the 
geodesies and therefore with the probability of presence of the 'particle'. Moreover, 



Q = -2V'—^ (7i 
JP 



can be interpreted as the new potential which is expected to emerge from the fractal 
geometry, in analogy with the identification of the gravitational field as a manifestation of 
the curved geometry in Einstein's general relativity. This result is supported by numerical 
simulations, in which the probability density is obtained directly from the distribution of 
geodesies without writing the Schrodinger equation [47, 103]. 



2.3.11 Nondifferentiable wave function 

In more recent works, instead of taking only the differentiable part of the velocity field 
into account, one constructs the covariant derivative and the wave function in terms of 
the full velocity field, including its divergent nondifferentiable part of zero mean [80, 94]. 
This still leads to the standard form of the Schrodinger equation. This means that, in 
the scale relativity framework, one expects the Schrodinger equation to have fractal and 
nondifferentiable solutions. This result agrees with a similar conclusion by Berry [10] and 
Hall [46], but it is considered here as a direct manifestation of the nondifferentiability of 
space itself. The research of such a behavior in laboratory experiments is an interesting 
new challenge for quantum physics. 
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2.4 Generalizations 



2.4.1 Fractal space time and relativistic quantum mechanics 

All these results can be generalized to relativistic quantum mechanics, that corresponds 
in the scale relativity framework to a full fractal space-time. This yields, as a first step, 
the Klein-Gordon equation [70, 72, 17]. 

Then the account of a new two-valuedness of the velocity allows one to suggest a 
geometric origin for the spin and to obtain the Dirac equation [17]. Indeed, the total 
derivative of a physical quantity also involves partial derivatives with respect to the space 
variables, d/dx^. From the very definition of derivatives, the discrete symmetry under 
the reflection dx^ ^ —dx^ is also broken. Since, at this level of description, one should 
also account for parity as in the standard quantum theory, this leads to introduce a bi- 
quaternionic velocity field [17], in terms of which Dirac bispinor wave function can be 
constructed. 

We refer the interested reader to the detailed papers [72, 17, 18]. 

2.4.2 Gauge fields as manifestations of fractal geometry 

The scale relativity principles has been also applied to the foundation of gauge theories, 
in the Abelian [70, 72] and non-Abelian [92, 103] cases. 

This application is based on a general description of the internal fractal structures of 
the "particle" (identified with the geodesies of a nondifferentiable space-time) in terms of 
scale variables rjap{x,y, z,t) = Qap^a^p whose true nature is tensorial, since it involves 
resolutions that may be different for the four space-time coordinates and may be corre- 
lated. This resolution tensor (similar to a covariance error matrix) generalizes the single 
resolution variable e. Moreover, one considers here a more profound level of description 
in which the scale variables may now be function of the coordinates. Namely, the internal 
structures of the geodesies may vary from place to place and during the time evolution, 
in agreement with the non-absolute character of the scale space. 

This generalization amounts to construct a 'general scale relativity' theory. The 
various ingredients of Yang-Mills theories (gauge covariant derivative, gauge invariance, 
charges, potentials, fields, etc..) can be recovered in such a framework, but they are now 
founded from first principles and are found to be of geometric origin, namely, gauge fields 
are understood as manifestations of the fractality of space-time [70, 72, 92, 103]. 

2.4.3 Quantum mechanics in scale space 

One may go still one step further, and also give up the hypothesis of differentiability of 
the scale variables. Another generalization of the theory then amounts to use in scale 
space the method that has been built for dealing with nondifferentiability in space-time 
[90] . This results in scale laws that take quantum-like forms instead of classical ones, and 
which may have several applications, as well in particle physics [90] as in biology [101]. 
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3 Applications 



3.1 Applications to physics and cosmology 

3.1.1 Application of special scale relativity: value of QCD coupling 

In the special scale relativity framework, the new status of the Planck length-scale as 
a lowest unpassable scale must be universal. In particular, it applies also to the de 
Broglie and Compton relations themselves. They must therefore be generalized, since in 
their standard definition they may reach the zero length, which is forbidden in the new 
framework. 

A fundamental consequence of these new relations for high energy physics is that the 
mass-energy scale and the length-time scale are no longer inverse as in standard quan- 
tum field theories, but they are now related by the special scale-relativistic generalized 
Compton formula, that reads [68] 

In ^ = (79) 
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As a consequence of this new relation, one finds that the grand unification scale be- 
comes the Planck energy scale [68, 69]. We have made the conjecture [68, 69] that the 
SU(3) inverse coupling reaches the critical value Att"^ at this unification scale, i.e., at an 
energy mpc^/27r in the special scale-relativistic modified standard model. 

By running the coupling from the Planck to the Z scale, this conjecture allows one 
to get a theoretical estimate for the value of the QCD coupling at Z scale. Indeed its 
renormalization group equation yields a variation of = as with length scale given to 
second order (for six quarks and Nh Higgs doublets) by 
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The variation with energy scale is obtained by making the transformation given by 
Eq. (79). This led in 1992 to the expectation [68] as^niz) = 0.1165±0.0005, that compared 
well with the experimental value at that time, a^lmz) = 0.112 ± 0.010, and was more 
precise. 

This calculation has been more recently reconsidered [12, 103], by using improved 
experimental values of the ai and a2 couplings at Z scale (which intervene at second 
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order), and by a better account of the top quark contribution. Indeed, its mass was 
unknown at the time of our first attempt in 1992, so that the running from Z scale to 
Planck scale was performed by assuming the contribution of six quarks on the whole scale 
range. 

However, the now known mass of the top quark, nit = 174.2 ± 3.3 GeV [108] is larger 
than the Z mass, so that only five quarks contribute to the running of the QCD coupling 
between Z scale and top quark scale, then six quarks between top and Planck scale. 
Moreover, the possibility of a threshold effect at top scale cannot be excluded. This led 
to an improved estimate : 

a,(mz) = 0.1173 ±0.0004, (81) 

which agrees within uncertainties with our initial estimate 0.1165(5) [68]. This expectation 
is in very good agreement with the recent experimental average as{mz) = 0.1176±0.0009 
[108], where the quoted uncertainty is the error on the average. We give in Fig. 1 the 
evolution of the measurement results of the strong coupling at Z scale, which compare 
very well with the theoretrical expectation. 
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Figure 1: Measured values of as{Mz) from 1992 (date of the theoretical prediction) to 2006 
[108] compared with the expectation as{mz) = 0.1173 it 0.0004 made from assuming that the 
inverse running coupling reaches the value 47r^ at Planck scale (see text). 
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3.1.2 Value of the cosmological constant 

One of the most difficuh open questions in present cosmology is the problem of the 
vacuum energy density and its manifestation as an effective cosmological constant. In 
the framework of the theory of scale relativity a new solution can be suggested to this 
problem, which also allows one to connect it to Dirac's large number hypothesis [69, Chap. 
7.1], [72]. 

The first step toward a solution has consisted in considering the vacuum as fractal, 
(i.e., explicitly scale dependent). As a consequence, the Planck value of the vacuum energy 
density is relevant only at the Planck scale, and becomes irrelevant at the cosmological 
scale. One expects such a scale-dependent vacuum energy density to be solution of a scale 
differential equation that reads 

dg/dlnr = r{g) = a + bg + 0{g^), (82) 

where g has been normalized to its Planck value, so that it is always < 1, allowing a 
Taylor expansion of T{g). This equation is solved as: 



Q = Qc 



l+(- 

r 



(83) 



This solution is the sum of a fractal, power law behavior at small scales, that can be 
identified with the quantum scale-dependent contribution, and of a scale- independent term 
at large scale, that can be identified with the geometric cosmological constant observed 
at cosmological scales. The new ingredient here is a fractal /non-fractal transition about 
some scale tq that comes out as an integration constant, and which allows to connect the 
two contributions. 

The second step toward a solution has been to realize that, when considering the 
various field contributions to the vacuum density, we may always chose < E >= (i.e., 
renormalize the energy density of the vacuum). But consider now the gravitational self- 
energy of vacuum fluctuations. It writes: 

^ G<E^> 

The Heisenberg relations prevent from making < E"^ >= 0, so that this gravitational 
self-energy cannot vanish. With < E"^ >i/2= hc/r, we obtain the asymptotic high energy 
behavior: 

Ir 



r 



^5) 



where is the Planck energy density and /p the Planck length. From this equation one 
can make the identification —b = 6, so that one obtains g = gc + (^o/''')^]- 
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Therefore one of Dirac's large number relations is proved from this result [69]. Indeed, 
introducing the characteristic length scale L = A~^/^ of the cosmological constant A 
(which is a curvature, i.e. the inverse of the square of a length), one obtains the relation: 

K = L//p = (ro//p)' = (mp/mo)^ (86) 

where the transition scale tq can be identified with the Compton length of a particle of 
mass mo. Then the power 3 in Dirac's large number relation is understood as coming 
from the power 6 of the gravitational self-energy of vacuum fluctuations and of the power 
2 that relies the invariant scale L to the cosmological constant, following the relation 
A = The important point here is that in this new form of the Eddington-Dirac's 

relation, the cosmological length is no longer the time- varying c/Hq (which led to theories 
of variation of constants), but the invariant cosmological length L, which can therefore 
be connected to an invariant elementary particle scale without no longer any need for 
fundamental constant variation. 

Now, a complete solution to the problem can be reached only provided the transition 
scale tq be identified. Our first suggestion [69, Chap. 7.1] has been that this scale is given 
by the classical radius of the electron. 

Let us give an argument in favor of this conjecture coming from a description of the 
evolution of the primeval universe. The classical radius of the electron actually defines 
the e~^e~ annihilation cross section and the e~e~ cross section a = irr"^ at energy mgC^. 
This length corresponds to an energy = Hc/tc = 70.02 MeV. This means that it yields 
the 'size' of an electron viewed by another electron. Therefore, when two electrons are 
separated by a distance smaller than re, they can no longer be considered as different, 
independent objects. 

The consequence of this property for the primeval universe is that should be a 
fundamental transition scale. When the Universe scale factor was so small that the inter- 
distance between any couple of electrons was smaller than re, there was no existing genuine 
separated electron. Then, when the cooling and expansion of the Universe separates the 
electron by distances larger than r^ the electrons that will later combine with the protons 
and form atoms appear for the first time as individual entities. Therefore the scale re and 
its corresponding energy 70 MeV defines a fundamental phase transition for the universe, 
which is the first appearance of electrons as we know them at large scales. Moreover, 
this is also the scale of maximal separation of quarks (in the pion), which means that the 
expansion, at the epoch this energy is reached, stops to apply to individual quarks and 
begins to apply to hadrons. This scale therefore becomes a reference static scale to which 
larger variable scales driven with the expansion can now be compared. Under this view, 
the cosmological constant would be a 'fossil' of this phase transition, in similarity with 
the 3K microwave radiation being a fossil of the combination of electrons and nucleons 
into atoms. 

One obtains with the CODATA 2002 values of the fundamental constants a theoretical 
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estimate 

K = (5.3000 ± 0.0012) x 10*^°, (87) 

i.e. Cu = InK = 139.82281(22), which corresponds to a cosmological constant (see [69] 
p. 305) 

A = (1.3628 ± 0.0004) x 10"^^ cm'^ (88) 
i.e., a scaled cosmological constant 

= (0.38874 ± 0.00012) h'^. (89) 
Finally the corresponding invariant cosmic length scale is theoretically predicted to be 

L = (2.77608 ± 0.00042) Gpc, (90) 

i.e., L = (8.5661 ± 0.0013) x 10^^ m. 

Let us compare these values with the most recent determinations of the cosmological 
constant (sometimes now termed, in a somewhat misleading way, 'dark energy'). The 
WMAP three year analysis of 2006 [123] has given h = 0.73 ± 0.03 and l]A(obs) = 
0.72 ± 0.03. These results, combined with the recent Sloan (SDSS) data [124], yield, 
assuming Q^ot = 1 (as supported by its WMAP determination, f2tot = 1-003 ± 0.010) 

Ac^ 

l]^(obs) = —2 = 0.761 ± 0.017, /i = 0.730 ±0.019. (91) 

3 -Ho 

Note that these recent results have also reinforced the cosmological constant interpretation 
of the 'dark energy' with a measurement of the coefficient of the equation of state w = 
—0.941 ± 0.094 [124], which encloses the value w = —1 expected for a cosmological 
constant. 

With these values one finds a still improved cosmological constant 

fiA/i2(obs) = 0.406 ± 0.030, (92) 

which corresponds to a cosmic scale 

L(obs) = (2.72 ± 0.10) Gpc, i.e., K(obs) = (5.19 ± 0.19) x lO*^", (93) 

in excellent agreement with the predicted values L(pred) = 2.7761(4) Gpc, and K = 
5.300(1) X 10^°. 

The evolution of these experimental determinations [103] is shown in Fig. 2 where they 
are compared with the theoretical expectation 

nAh^ipred) = 0.38874 ± 0.00012. (94) 

The convergence of the observational values toward the theoretical estimate, despite an 
improvement of the precision by a factor of more than 20, is striking. The 2008 value 
from the Five- Year WMAP results is r2A^^(obs) = 0.384 ± 0.043 [49] and is once again in 
very good agreement with the theoretical expectation made 16 years ago [69], before the 
first genuine measurements in 1998. 
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Figure 2: Evolution of the measured values of the dimensionless cosmological constant Q.\h'^ = 
Ac^/3i/^QQ, from 1975 to 2006, compared to the theoretical expectation A = (me/amp)^ (1/^p)^ 
[69] that gives numerically OA^^(pred) = 0.38874 ± 0.00012. 



3.2 Applications to astrophysics 
3.2.1 Gravitational Schrodinger equation 

Let us first briefly recall the basics of the scale-relativistic theoretical approach. It has 
been reviewed in Sec. 2.3 in the context of the foundation of microphysics quantum me- 
chanics. We shall now see that some of its ingredients, leading in particular to obtain 
a generalized Schrodinger form for the equation of motion, also applies to gravitational 
structure formation. 

Under three general conditions, namely, {(i) infinity of geodesies (which leads to in- 
troduce a non-deterministic velocity field), (ii) fractal dimension Dp = 2 oi each geodesic, 
on which the elementary displacements are described in terms of the sum dX = dx + 
of a classical, differentiable part dx and of a fractal, non-differentiable fluctuation c?^, (iii) 
two-valuedness of the velocity field, which is a consequence of time irreversibility at the 
infinitesimal level issued from non-differentiability, one can construct a complex covariant 
derivative that reads ^ 

| = | + V.V->I,A, (95) 
where D is a parameter that characterizes the fractal fluctuation, which is such that 
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< d^'^ >= 2Vdt, and where the classical part of the velocity field, V is complex as a 
consequence of condition (iii) (see [17, 96] for more complete demonstrations). 

Then this covariant derivative, that describes the non-differentiable and fractal geome- 
try of space-time, can be combined with the covariant derivative of general relativity, that 
describes the curved geometry. We shall briefiy consider in what follows only the Newto- 
nian limit. In this case the equation of geodesies keeps the form of Newton's fundamental 
equation of dynamics in a gravitational field, 

SV SV vflUo. (96) 



dt dt KfTT', 

where (j) is the Newtonian potential energy. Introducing the action S, which is now 
complex, and making the change of variable ^ = e*'^/^'"-'^, this equation can be integrated 
under the form of a generalized Schrodinger equation [69]: 

V'^AiP + iV^iP - —iP = 0. (97) 
ot 2m 

Since the imaginary part of this equation is the equation of continuity (Sec. 3), and 
basing ourselves on our description of the motion in terms of an infinite family of geodesies, 
P = naturally gives the probability density of the particle position [17, 96]. 

Even though it takes this Schrodinger-like form, equation (97) is still in essence an 
equation of gravitation, so that it must come under the equivalence principle [73, 2], 
i.e., it is independent of the mass of the test-particle. In the Kepler central potential case 
(0 = —GMm/r), GM provides the natural length-unit of the system under consideration. 
As a consequence, the parameter V reads: 

GM , , 

where w is a constant that has the dimension of a velocity. The ratio Ug = w/c actually 
plays the role of a macroscopic gravitational coupling constant [2, 81]. 



3.2.2 Formation and evolution of structures 

Let us now compare our approach with the standard theory of gravitational structure 
formation and evolution. By separating the real and imaginary parts of the Schrodinger 
equation we obtain, after a new change of variables, respectively a generalized Euler- 
Newton equation and a continuity equation, namely, 

m{^^+V-V)V = -V(0 + g), ^ + div(PV) = 0, (99) 

where V is the real part of the complex velocity field V. In the case when the density of 
probability is proportional to the density of matter, P oc p, this system of equations is 
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equivalent to the classical one used in the standard approach of gravitational structure 
formation, except for the appearance of an extra potential energy term Q that writes: 



The existence of this potential energy, (which amount to the Bohm potential in stan- 
dard quantum mechanics) is, in our approach, readily demonstrated and understood: 
namely, it is the very manifestation of the fractality of space, in similarity with New- 
ton's potential being a manifestation of curvature. We have suggested [82, 91, 93] that it 
could be the origin of the various effects which are usually attributed to an unseen, 'dark' 
matter. 

In the case when actual particles achieve the probability density distribution (structure 
formation), we have p = moP; then the Poisson equation (i.e., the field equation) becomes 
A0 = ATrGmmolipl'^ and it is therefore strongly interconnected with the Schrodinger equa- 
tion (which is here a new form for the equation of motion). Such a system of equations 
is similar to that encountered in the description of superconductivity (Hartree equation). 
We expect its solutions to provide us with general theoretical predictions for the structures 
(in position and velocity space) of self-gravitating systems at multiple scales [74, 27]. This 
expectation is already supported by the observed agreement of several of these solutions 
with astrophysical observational data [69, 73, 81, 76, 79, 77, 78, 48]. 

3.2.3 Planetary systems 

Let us briefly consider the application of the theory to the formation of planetary systems. 
The standard model of formation of planetary systems can be reconsidered in terms of 
a fractal description of the motion of planetesimals in the protoplanetary nebula. On 
length-scales much larger than their mean free path, we have assumed [69] that their highly 
chaotic motion satisfy the three conditions upon which the derivation of a Schrodinger 
equation is based (large number of trajectories, fractality and time symmetry breaking). 
In modern terms, our proposal is but a 'migration' theory, since it amounts to take into 
account the coupling between planetesimals (or proto-planets) and the remaining disk. 
But, instead of considering a mean field coupling, we consider the effect of the closest 
bodies to be the main one, leading to Brownian motion and irreversibility. 

This description applies to the distribution of planetesimals in the proto-planetary 
nebula at several embedded levels of hierarchy. Each hierarchical level (k) is characterized 
by a length-scale defining the parameter Vk (and therefore the velocity Wk) that appears in 
the generalized Schrodinger equation describing this sub-system. This hierarchical model 
has allowed us to recover the mass distribution of planets and small planets in the inner 
and outer solar systems [76]. It is generally supported by the structure of our own solar 
system, which is made of several subsystems embedded one in another, namely: 



Q 




(100) 
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The Sun Through Kepler's third law, the velocity w = 3 x 144.7 = 434.1 km/s is 
very closely the Keplerian velocity at the Sun radius {Rq = 0.00465 AU corresponds to 
Wq = 437.1 km/s). Moreover, one can also apply our approach to the organization of the 
sun surface itself. One expect the distribution of the various relevant physical quantities 
that characterize the solar activity at the Sun surface (sun spot number, magnetic field, 
etc..) to be described by a wave function whose stationary solutions read ip = ipo e*-^*/^™-^. 

The energy E results from the rotational velocity and, to be complete, should also 
include the turbulent velocity, so that E = {v'^^^ + v^^j.^^)/2. This means that we expect 
the solar surface activity to be subjected to a fundamental period: 

27rmV A-kV 

^rot + ^'turb 

The parameter T) at the Sun radius is P = GMq/2wq, then we obtain: 



'r = —F^ = -2 (101) 



^ = 102 

The average sideral rotation period of the Sun is 25.38 days, yielding a velocity of 2.01 
km/s at equator [109]. The turbulent velocity has been found to be fturb = 1.4±0.2 km/s 
[56]. Therefore we find numerically 

r = (10.2 ± 1.0) yrs. (103) 

The observed value of the period of the Solar activity cycle, robs = 11.0 yrs, nicely supports 
this theoretical prediction. This is an interesting result, owing to the fact that there is, 
up to now, no existing theoretical prediction of the value of the solar cycle period, except 
in terms of very rough order of magnitude [128]. 

Moreover, since we have now at our disposal a simple and precise formula for a stellar 
cycle which precisely accounts for the solar period, the advantage is that it can be tested 
with other stars. The observation of the magnetic activity cycle of distant solar-like stars 
remains a difficult task, but it has now been performed on several stars. A first attempt 
gives very encouraging results (see Fig. 3), since we obtain indeed a satisfactory agreement 
between the observed and predicted periods, in a statistically significant way, despite the 
small number of objects. 



The intramercurial system organized on the constant Wq = 3 x 144 = 432 km/s. 
The existence of an intramercurial subsystem is supported by various stable and transient 
structures observed in dust, asteroid and comet distributions (see [27]). We have in 
particular suggested the existence of a new ring of asteroids, the 'Vulcanoid belt', at a 
preferential distance of about 0.17 AU from the Sun. 



The inner solar system (earth-like planets), organized with a constant Wi = 144 km/s 
(see Fig. 6). 
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Figure 3: Comparison between the observed values of the period of solar-hke star cycles (inactive 
stars with better determined behavior in Table 1 of Ref. [116]) and the predicted periods (see 
text). The open point is for the Sun. The correlation is significant at a probability level P ~ 10~^ 
(Student variable t 5). 



The outer solar system organized with a constant Wo = 144/5 = 29 km/s (see Fig. 4), 
as deduced from the fact that the mass peak of the inner solar system lies at the Earth 
distance (n = 5). The Jovian giant planets lie from n = 2 to n = 5. Pluton lies on n = 6, 
but is now considered to be a dwarf planet part of the Kuiper belt. 

Kuiper belt The recently discovered Kuiper and scattered Kuiper belt objects (Fig. 
4) show peaks of probability at = 6 to 9 [27], as predicted before their discovery [71]. 
In particular, the predicted peak around 57 AU {n = 7) is the main observed peak in 
the SKBO distribution of semi- major axes. The following peak [n = 8), predicted to 
be around 70 AU, has received a spectacular verification with the discovery of the dwarf 
planet Eris (2003 UB313) at 68 AU, whose mass larger than Pluton has recently led to a 
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Figure 4: Distribution of the semi-major axis of Kuiper belt objects (KBO) and scattered 
Kuiper belt objects (SKBO), compared with the theoretical predictions (arrows) of probability 
density peaks for the outer solar system [27] (see text). The existence of probability density 
peaks for the Kuiper belt at f» 40, 55, 70, 90 AU, etc., has been theoretically predicted in 1993 
before the discovery of these objects [71], and it is now supported by the observational data, in 
particular by the new small planet Eris at 68 AU, whose mass is larger than Pluto, and which 
falls close to the expected probability peak n = 8 at 70 AU (see text). 



revision of planetary nomenclature. 

Distant Kuiper belt Beyond these distances, we have been able to predict a new 
level of hierarchy in the Solar System whose main SKBO peak at 57 AU would be the 
fundamental level {n = 1) [37]. The following probability peaks are expected, accord- 
ing to the law, to lie for semi-major axes of 228, 513, 912, 1425, 2052 AU, etc.... 
Once again this prediction has been validated by the observational data in a remark- 
able way (see Fig. 5), since 4 bodies, including the very distant small planet Sedna, 
have now been discovered in the 513 AU peak {n = 3), 7 bodies in the 228 AU peak 
{n = 2) , and now one very distant object at about 1000 AU (data Minor Planet Center, 
http:/ /www. cfa.harvard.edu/iau/mpc. html). 
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Figure 5: Distribution of the semi-major axis of very distant scattered Kuiper belt objects 
(SKBO) , compared with the theoretical predictions of probability density peaks (see text). We 
have taken the main SKBO peak at ~57 AU (which is the predicted n = 7 peak of the outer 
solar system) as fundamental level (n = 1) for this new level of hierarchy of the solar system. 
The figure plots the histogram of the variable (o/57)^/^, where a is the semimajor axis of the 
object orbit in AU. The theoretical prediction, done before the discovery of the distant objects, 
is that the distribution of this variable should show peaks for integer values, as now verified by 
the observational data. 



Extrasolar planets We have suggested more than 16 years ago [69, 71], before the 
discovery of exoplanets, that the theoretical predictions from this approach of planetary 
formation should apply to all planetary systems, not only our own solar system. Mean- 
while more than 300 exoplanets have now been discovered, and the observational data 
support this prediction in a highly statistically significant way (see [73, 81, 27] and Fig. 6). 

The presently known exoplanets mainly correspond to the intramercurial and inner 
solar systems. The theoretical prediction, made in 1993 [69, Chap. 7.2], according to 
which the distribution of semi-major axes a is expected to show peaks of probability for 
integer values of the variable 4.83 (a/M)^/^, where M is the star mass, remains validated 
with a high statistical significance (see Fig. 6). In particular, in addition to the peaks 
of probability corresponding to the inner solar system planets (n = 3 Mercury, n = 4 
Venus, n = 5 Earth, n = Q Mars), two additional predicted peaks of probability, the 
'fundamental' one at 0.043 AU/M© and the second one at 0.17 AU/Mq, have been made 
manifest in extrasolar planetary systems. In particular, the validation of the principal 
prediction of the SR approach, namely, the main peak at the fundamental level ra = 1, is 
striking since it now contains more than 80 exoplanets. 
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Figure 6: Observed distribution of the semi-major axes of 300 exoplanets (June 2008 data [118]) 
and inner solar system planets, compared with the theoretical prediction. The figure gives the 
histogram of the distribution of the variable 4.83(a/M)^/^. One predicts the occurence of peaks 
of probabihty density for semimajor axes a„ = GM^u/wq)^, where n is integer, M is the star 
mass and wq = 144.7 it 0.7 km/s is a gravitational coupling constant (see text). The main peak 
at n = 1 (fundamental level), which now contains 74 exoplanets, has been cut to better view the 
secondary peaks. The probability to obtain such an agreement by chance between the predicted 
(red vertical lines) and observed peaks is now P = 5 x 10~^. 



3.3 Applications to sciences of life 

The scale relativity theory has also been recently applied to sciences other than physical 
sciences, including sciences of life, sciences of societies, historical [34] and geographical 
sciences [61, 35, 36] and human sciences [125, 87, 95, 102]. We refer the interested reader 
to the books [85, 99], to parts of review papers or books [83, 90, 97] and full review papers 
on this specific subject [6, 101] for more details. 

3.3.1 Applications of log-periodic laws 

Species evolution Let us first consider the application of log-periodic laws to the 
description of critical time evolution. Recall that a log-periodic generalization to scale 
invariance has been obtained as a solution to wave-like differential scale equations. Inter- 
preted as a distribution of probability, such solutions therefore lead to a geometric law of 
progression of probability peaks for the occurence of events. 

Now several studies have shown that many biological, natural, sociological and eco- 
nomic phenomena obey a log-periodic law of time evolution such as can be found in some 
critical phenomena : earthquakes [120], stock market crashes [121], evolutionary leaps 
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Figure 7: The dates of major evolutionary events of seven lineages (common evolution from life 
origin to viviparity, Theropod and Sauropod dinosaurs, Rodents, Equidae, Primates including 
Hominidae, and Echinoderms) are plotted as black points in terms of log(Tc — T), and compared 
with the numerical values from their corresponding log-periodic models (computed with their 
best-fit parameters). The adjusted critical time Tc and scale ratio g are indicated for each lineage 
[19, 85, 86]. 



[19, 85, 86], long time scale evolution of western and other civilizations [85, 86, 45], world 
economy indices dynamics [50], embryogenesis [14], etc... Thus emerges the idea that 
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this behaviour typical of temporal crisis could be extremely widespread, as much in the 
organic world as in the inorganic one [122]. 

In the case of species evolution, one observes the occurrence of major evolutionary 
leaps leading to bifurcations among species, which proves the existence of punctuated 
evolution [41] in addition to the gradual one. The global pattern is assimilated to a 'tree 
of life', whose bifurcations are identified to evolutionary leaps, and branch lengths to 
the time intervals between these major events [19]. As early recognized by Leonardo da 
Vinci, the branching of vegetal trees and rivers may be described as a first self-similar 
approximation by simply writing that the ratio of the lengths of two adjacent levels is 
constant in the mean. We have made a similar hypothesis for the time intervals between 
evolutionary leaps, namely, (T„ — T„_i)/(T„_|_i — T„) = g. Such a geometric progression 
yields a log-periodic acceleration for g > 1, a, deceleration for g <!■, and a periodicity for 
g = 1. Except when g = 1, the events converge toward a critical time which can then 
be taken as reference, yielding the following law for the event T„ in terms of the rank n: 



where Tq is any event in the lineage, n the rank of occurrence of a given event and g is 
the scale ratio between successive time intervals. Such a chronology is periodic in terms 
of logarithmic variables, i.e., log |T„ — = log \Tq — — nlogg. 

This law is dependent on two parameters only, g and Tc, which of course have no 
reason a priori to be constant for the entire tree of life. Note that g is not expected to 
be an absolute parameter, since it depends on the density of events chosen, i.e., on the 
adopted threshhold in the choice of their importance (namely, if the number of events 
is doubled, g is replaced by ^/g)■ Only a maximal value of g, corresponding to the very 
major events, could possibly have a meaning. On the contrary, the value of Tc is expected 
to be a characteristic of a given lineage, and therefore not to depend on such a choice. 
This expectation is supported by an analysis of the fossil record data under various choices 
of the threshold on the events, which have yielded the same values of within error bars 



A statistically significant log-periodic acceleration has been found at various scales 
for global life evolution, for primates, for sauropod and theropod dinosaurs, for rodents 
and North American equids. A deceleration law was conversely found in a statistically 
significant way for echinoderms and for the first steps of rodents evolution (see Fig. 7 and 
more detail in Refs. [19, 85, 86]). One finds either an acceleration toward a critical date 
Tc or a deceleration from a critical date, depending on the considered lineage. 

It must be remarked that the observed dates follow a log-periodic law only in the 
mean, and show a dispersion around this mean (see [85, p. 320]. In other words, this 
is a statistical acceleration or deceleration, so that the most plausible interpretation is 
that the discrete T„ values are nothing but the dates of peaks in a continuous probability 
distribution of the events. Moreover, it must also be emphasized that this result does not 
put the average constancy of the mutation rate in question. This is demonstrated by a 
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[86]. 
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study of the cytochrome c tree of branching (in preparation), which is based on genetic 
distances instead of geological chronology, and which nevertheless yields the same result, 
namely, a log-periodic acceleration of most lineages, and a periodicity (which corresponds 
to a critical time tending to infinity) in some cases. The average mutation rate remains 
around 1/20 Myr since about 1 Gyr, so that one cannot escape the conclusion that the 
number of mutations needed to obtain a major evolutionary leap decreases with time 
among many lineages, and increases for some of them. 

Embryogenesis and human development Considering the relationships between 
phylogeny and ontogeny, it appeared interesting to verify whether the log-periodic law 
describing the chronology of several lineages of species evolution may also be applied 
to the various stages in human embryological development. The result, (see Figure in 
Chaline's contribution), is that a statistically significant log-periodic deceleration with a 
scale ratio g = 1.71 ±0.01 is indeed observed, starting from a critical date that is consitent 
with the conception date [14]. 

Evolution of societies Many observers have commented on the way historical events 
accelerate. Grou [44] has shown that the economic evolution since the neolithic can be 
described in terms of various dominating poles which are submitted to an accelerating 
crisis-nocrisis pattern, which has subsequently been quantitatively analysed using log- 
periodic laws. 

For the Western civilization since the Neolithic (i.e., on a time scale of about 8000 
years), one finds that a log-periodic acceleration with scale factor g = 1.32 ±0.018 occurs 
toward Tc = 2080 ±30 (see Fig. 8), in a statistically highly significant way. This result has 
been later confirmed by Johansen and Sornette [50] by an independent study on various 
market, domestic, research and development, etc... indices on a time scale of about 200 
years, completed by demography on a time scale of about 2000 years. They find critical 
dates for these various indices in the range 2050-2070, which support the longer time scale 
result. 

One of the intriguing features of all these results is the frequent occurence of values 
of the scale ratio g close to g = 1.73 and its square root 1.32 (recall that one passes from 
a value of g to its square root by simply doubling the number of events). This suggests 
once again a discretization of the values of this scale ratio, that may be the result of a 
probability law (in scale space) showing quantized probability peaks . We have considered 
the possibility that (7 = 1.73 ~ \/3 could be linked to a most probable branching ratio of 3 
[19, 85], while Queiros-Conde [113] has proposed a 'fractal skin' model for understanding 
this value. 
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Figure 8: Comparison of the median dates of the main economic crises of western civiHzation 
with a log-periodic accelerating law of critical date Tc = 2080 and scale ratio g = 1.32 (figm'e 
a). The last white point corresponds to the predicted next crisis (1997-2000) at the date of the 
study (1996), as has been later supported in particular by the 1998 and 2000 market crashes, 
while the next crises are now predicted for (2015-2020), then (2030-2035). Figure b shows the 
estimation of the critical date through the optimisation of the Student's t variable. This result is 
statistically significant, since the probability to obtain such a high peak by chance is P < 10~^. 



3.3.2 History and geography 

The application of the various tools and methods of the scale relativity theory to history 
and geography has been proposed by Martin and Forriez [34, 61, 35, 36]. Forriez has 
shown that the chronology of some historical events (various steps of evolution of a given 
site) recovered from archeological and historical studies can be fitted by a log-periodic 
deceleration law with again g ^ 1.7 and a retroprediction of the foundation date of the site 
from the critical date [34, 35]. Moreover, the various differential equation tools developed 
in the scale relativity approach both in scale and position space, including the nonlinear 
cases of variable fractal dimensions, have been found to be particularly well adapted to 
the solution of geographical problems [61]. 

3.3.3 Predictivity 

Although these studies remain, at that stage, of an empirical nature (it is only a purely 
chronological analysis which does not take into account the nature of the events), they 
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nevertheless provide us with a beginning of predictivity. Indeed, the fitting law is a two 
parameter function (Tc and g) that is applied to time intervals, so that only three events 
are needed to define these parameters. Therefore the subsequent dates are predicted after 
the third one, in a statistical way. Namely, as already remarked, the predicted dates 
should be interpreted as the dates of the peaks of probability for an event to happen. 
Examples of such a predictivity (or retropredictivity) are: 

(i) the retroprediction that the common Homo- Pan-Gorilla ancestor (expected, e.g., from 
genetic distances and phylogenetic studies), has a more probable date of appearance at 
^ —10 millions years [19]; its fossil has not yet been discovered (this is one of the few 
remaining 'missing links'); 

(ii) the prediction of a critical date for the long term evolution of human societies around 
the years 2050-2080 [85, 50, 86, 45]; 

(iii) the finding that the critical dates of rodents may reach +60 Myrs in the future, 
showing their large capacity of evolution, in agreement with their known high biodiversity; 

(iv) the finding that the critical dates of dinosaurs are about —150 Myrs in the past, 
indicating that they had reached the end of their capacity of evolution (at least for the 
specific morphological characters studied) well before their extinction at —65 Myrs; 

(v) the finding that the critical dates of North american Equids is, within uncertain- 
ties, consistent with the date of their extinction, which may mean that, contrarily to the 
dinosaur case, the end of their capacity of evolution has occured during a phase of en- 
vironmental change that they have not been able to deal with by the mutation-selection 
process; 

(vi) the finding that the critical date of echinoderms (which decelerate instead of accelerat- 
ing) is, within uncertainties, the same as that of their apparition during the PreCambrian- 
Cambrian radiation, this supporting the view of the subsequent events as a kind of "scale 
wave" expanding from this first shock. 

3.3.4 Applications in Earth sciences 

As last examples of such a predictivity, let us give some examples of applications of 
critical laws (power laws in |T — T^p and their log-perodic generalizations) to problems 
encountered in Earth sciences, namely, earthquakes (California and Sichuan) and decline 
of Arctic sea ice. 

California earthquakes The study of earthquakes has been one of the first domain of 
application of critical and log-periodic laws [120, 3]. The rate of California earthquakes 
is found to show a very marked log-periodic deceleration [98, 99]. We show indeed in 
Fig. 9 the observed rate of Southern California earhquakes of magnitude larger than 5, 
compared with a log-periodic deceleration law. This model allows us to predict future 
peaks of probability around the years 2050 then 2115. 



43 




1800 1850 1900 1950 2000 2050 

Date 

Figure 9: Observed rate of Southern California earhquakes of magnitude larger than 5 (his- 
togram). The data are taken from the U.S. Geological Survey EarthQuake Data Center (years 
1932-2006) and EarthQuake Data Base (Historical earthquakes, years 1500-1932). This rate is 
well fitted by a power law subjected to a log-periodic fluctuation decelerating since a critical 
date Tc = 1796 (red fluctuating line). The model predicts the next probability peak around the 
years 2050 [98]. 



Sichuan 2008 earthquake The May 2008 Sichuan earthquake and its replicas also 
yields a good example of log-periodic deceleration, but on a much smaller time scale (see 
Fig 10). 

Arctic sea ice extent It is now well-known that the decrease of arctic sea-ice extent has 
shown a strong acceleration in 2007 and 2008 with respect to the current models assuming 
a constant rate (~ 8% by decade), which predicted in 2006 a total disappearance of the 
ice at minimum (15 September) for the end of the century. From the view point of these 
models, the 2007 and now 2008 values (see Fig. 11) were totally unexpected. 

However, we have proposed, before the knowledge of the 2007 minimum, to fit the 
data with a critical law of the kind y = i/q — a\T — Td'^ . Such an accelerating law 
has the advantage to include in its structure the fact that one expect the ice to fully 
disappear after some date, while the constant rate law formally pushed the date of total 
disappearance to infinity. The fit of the data up to 2006 with the critical law was already 
far better than with the constant rate law, and it actually allowed us to predict a full 
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Figure 10: Time evolution during 14 days of the replicas of the May 12, 2008 Sichuan earth- 
quake (data obtained and studied May 27, 2008 from the seismic data bank EduSeis Explorer, 
http://aster.unice.fr/EduSeisExplorer/form-sis.asp). The (up) figure gives the magnitudes of 
the replicas and the (down) figure the rate of replicas. Both show a continuous decrease to 
which are added discrete sharp peaks. The peaks which are common to both diagrams show a 
clear deceleration according to a log-periodic law starting from the main earthquake (red ver- 
tical lines), which allows one to predict the next strongest replicas with a good precision. For 
example, the peak of replicas of 25 May 2008 could be predicted with a precision of 1.5 day from 
the previous peaks. Reversely, the date of the main earthquake (May 12.27 2008, magnitude 
7.9) can be retropredicted from that of the replicas with a precision of 6 h. 



disappearance epoch far closer than previously expected and a low 2007 point [100]. The 
2008 point has confirmed the validity of the model in an impressive way (Fig. 11). We 
obtain by the method a best fit for the minimum ice surface (in square kilometers), 
y = 8 — 12.3 X |T — 2012|~°'*^. The critical time is as early as Tc = 2012, which means 



45 



U 

P 2 - 
w 

1 - 



1980 1985 1990 1995 2000 2005 2010 

Date 

Figure 11: Observed evolution of the minimum arctic sea ice extent, according to the data of 
the U.S. National Snow and Ice Data Center (NSIDC, http://nsidc.org/), from 1979 to 2008. 
The minimum ocurs around 15 September of each year. This evolution is compared to: (i) the 
standard fit corresponding to an assumed constant rate of extent decrease (blue line) ; (ii) a fit 
by a critical law accelerating toward a critical date Tc = 2012. The second fit is far better and 
has allowed us to predict the 2007 and 2008 low points before their observation [100]. It implies 
that the arctic sea is expected to be totally free from ice by September 2011. 



that a full ice melting is predicted for September 2011, and is even possible for September 
2010, for which the model gives only 1.2 million km^ of remaining ice surface. 

The application of the same method to the mean surface data during August and 
October months also shows a clear acceleration toward Tc = 2013, which means that only 
one year later (2012) the arctic sea can be expected to be free from ice during several 
months (August to October). 

3.4 Applications of scale relativity to biology 

One may consider several applications to biology of the various tools and methods of the 
scale relativity theory, namely, generalized scale laws, macroscopic quantum-type theory 
and Schrodinger equation in position space then in scale space and emergence of gauge- 
type fields and their associated charges from fractal geometry [69, 85, 89, 6, 96]. One knows 
that biology is founded on biochemistry, which is itself based on thermodynamics, to which 
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we contemplate the future possibility to apply the macroquantization tools described in 
the theoretical part of this article. Another example of future possible applications is to 
the description of the growth of polymer chains, which could have consequences for our 
understanding of the nature of DNA and RNA molecules. 
Let us give some explicit examples of such applications. 

3.5 Confinement 

The solutions of non-linear scale equations such as that involving a harmonic oscillator- 
like scale force [75] may be meaningful for biological systems. Indeed, its main feature 
is its capacity to describe a system in which a clear separation has emerged between an 
inner and an outer region, which is one of the properties of the first prokaryotic cell. We 
have seen that the effect of a scale harmonic oscillator force results in a confinement of 
the large scale material in such a way that the small scales may remain unaffected. 

Another interpretation of this scale behavior amounts to identify the zone where the 
fractal dimension diverges (which corresponds to an increased 'thickness' of the material) 
as the description of a membrane. It is indeed the very nature of biological systems to have 
not only a well-defined size and a well-defined separation between interior and exterior, 
but also systematically an interface between them, such as membranes or walls. This 
is already true of the simplest prokaryote living cells. Therefore this result suggests the 
possibility that there could exist a connection between the existence of a scale field (e.g., 
a global pulsation of the system, etc..) both with the confinement of the cellular material 
and with the appearance of a limiting membrane or wall [89]. This is reminiscent of 
eukaryotic cellular division which involves both a dissolution of the nucleus membrane and 
a deconfinement of the nucleus material, transforming, before the division, an eukaryote 
into a prokaryote-like cell. This could be a key toward a better understanding of the first 
major evolutionary leap after the appearance of cells, namely the emergence of eukaryotes. 

3.6 Morphogenesis 

The generalized Schrodinger equation (in which the Planck constant h can be replaced 
by a macroscopic constant) can be viewed as a fundamental equation of morphogenesis. 
It has not been yet considered as such, because its unique domain of application was, up 
to now, the microscopic (molecular, atomic, nuclear and elementary particle) domain, in 
which the available information was mainly about energy and momentum. 

However, scale relativity extends the potential domain of application of Schrodinger- 
like equations to every systems in which the three conditions (infinite or very large number 
of trajectories, fractal dimension of individual trajectories, local irreversibility) are ful- 
filled. Macroscopic Schrodinger equations can be constructed, which are not based on 
Planck's constant h, but on constants that are specific of each system (and may emerge 
from their self-organization). 
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Now the three above conditions seems to be particularly well adapted to the description 
of living systems. Let us give a simple example of such an application. 

In living systems, morphologies are acquired through growth processes. One can 
attempt to describe such a growth in terms of an infinite family of virtual, fractal and 
locally irreversible, trajectories. Their equation can therefore be written under the form 
of a fractal geodesic equation, then it can be integrated as a Schrodinger equation. 

If one now looks for solutions describing a growth from a center, one finds that this 
problem is formally identical to the problem of the formation of planetary nebulae [27] , 
and, from the quantum point of view, to the problem of particle scattering, e.g., on an 
atom. The solutions looked for correspond to the case of the outgoing spherical probability 
wave. 

Depending on the potential, on the boundary conditions and on the symmetry con- 
ditions, a large family of solutions can be obtained. Considering here only the simplest 
ones, i.e., central potential and spherical symmetry, the probability density distribution of 
the various possible values of the angles are given in this case by the spherical harmonics. 



These functions show peaks of probability for some angles, depending on the quantized 
values of the square of angular momentum (measured by the quantum number /) and 
of its projection on axis z (measured by the quantum number m). 




Figure 12: Morphogenesis of a 'flower'-like structure, solution of a generalized Schrodinger 
equation that describes a growth process from a center(/ = 5, m = 0). The 'petals', 'sepals' and 
'stamen' are traced along angles of maximal probability density. A constant force of 'tension' 
has been added, involving an additional curvature of 'petals', and a quantization of the angle 6 
that gives an integer number of 'petals' (here, fc = 5). 

Finally a more probable morphology is obtained by 'sending' matter along angles of 
maximal probability. The biological constraints leads one to skip to cylindrical symmetry. 
This yields in the simplest case a periodic quantization of the angle 9 (measured by an 
additional quantum number /c), that gives rise to a separation of discretized 'petals'. 




(105) 
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Moreover there is a discrete symmetry breaking along the z axis hnked to orientation 
(separation of 'up' and 'down' due to gravity, growth from a stem). The solutions obtained 
in this way show floral 'tulip '-like shapes (see Fig. 12 and [83, 89, 96]). 

Coming back to the foundation of the theory, it is remarkable that these shapes are 
solutions of a geodesic, strongly covariant equation dV/dt = 0, which has the form of 
the Galilean motion equation in vacuum in the absence of external force. Even more 
profoundly, this equation does not describe the motion of a particle, but purely geometric 
virtual possible paths, this given rise to a description in terms of a probability density 
which plays the role of a potential for the real particle (if any, since, in the application to 
elementary particles, we identify the 'particles' with the geodesies themselves, i.e., they 
become pure relative geometric entities devoid of any proper existence) . 

3.7 Origin of life 

The problems of origin are in general more complex than the problems of evolution. 
Strictly, there is no 'origin' and both problems could appear to be similar, since the 
scientific and causal view amounts to consider that any given system finds its origin in an 
evolution process. However, systems are in general said to evolve if they keep their nature, 
while the question is posed in terms of origin when a given system appears from another 
system of a completely different nature, and moreover, often on times scales which are 
very short with respect to the evolution time. An example in astrophysics is the origin 
of stars and planetary systems from the interstellar medium, and in biology the probable 
origin of life from a prebiotic medium. 

A fondamentally new feature of the scale relativity approach concerning such problems 
is that the Schrodinger form taken by the geodesic equation can be interpreted as a 
general tendency for systems to which it applies to make structures, i.e., to lead to self- 
organization. In the framework of a classical deterministic approach, the question of 
the formation of a system is always posed in terms of initial conditions. In the new 
framework, the general existence of stationary solutions allows structures to be formed 
whatever the initial conditions, in correspondence with the field, the symmetries and the 
boundary conditions (namely the environmental conditions in biology), and in function 
of the values of the various conservative quantities that characterize the system. 

Such an approach could allow one to ask the question of the origin of life in a renewed 
way. This problem is the analog of the 'vacuum' (lowest energy) solutions, i.e., of the 
passage from a non-structured medium to the simplest, fundamental level structures. In 
astrophysics and cosmology, the problem amounts to understand the apparition, from the 
action of gravitation alone, of structures (planets, stars, galaxies, clusters of galaxies, large 
scale structures of the Universe) from a highly homogeneous and non-structured medium 
whose relative fluctuations were smaller than 10~^ at the time of atom formation. In the 
standard approach to this problem a large quantity of postulated and unobserved dark 
matter is needed to form structures, and even with this help the result is dissatisfying. In 
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the scale relativity framework, we have suggested that the fundamentally chaotic behavior 
of particle trajectories leads to an underlying fractal geometry of space, which involves 
a Schrodinger form for the equation of motion, leading both to a natural tendency to 
form structures and to the emergence of an additional potential energy, identified with 
the 'missing mass(-energy)'. 

The problem of the origin of life, although clearly far more difficult and complex, 
shows common features with this question. In both cases one needs to understand the 
apparition of new structures, functions, properties, etc... from a medium which does not 
yet show such structures and functions. In other words, one need a theory of emergence. 
We hope that scale relativity is a good candidate for such a theory, since it owns the two 
required properties: (i) for problems of origin, it gives the conditions under which a weakly 
structuring or destructuring (e.g., diffusive) classical system may become quantum-like 
and therefore structured; (ii) for problems of evolution, it makes use of the self-organizing 
property of the quantum-like theory. 

We therefore tentatively suggest a new way to tackle the question of the origin of life 
(and in parallel, of the present functionning of the intracellular medium) [89, 6, 101]. The 
prebiotic medium on the primordial Earth is expected to have become chaotic in such a 
way that, on time scales long with respect to the chaos time (horizon of predictibility) , the 
conditions that underlie the transformation of the motion equation into a Schrodinger- 
type equation, namely, complete information loss on angles, position and time leading 
to a fractal dimension 2 behavior on a range of scales reaching a ratio of at least 10^- 
10^, be fulfilled. Since the chemical structures of the prebiotic medium have their lowest 
scales at the atomic size, this means that, under such a scenario, one expects the first 
organized units to have appeared at a scale of about 10 /im, which is indeed a typical 
scale for the first observed prokaryotic cells. The spontaneous transformation of a classical, 
possibly diffusive mechanics, into a quantum-like mechanics, with the diffusion coefficient 
becoming the quantum self-organization parameter T) would have immediate dramatic 
consequences: quantization of energy and energy exchanges and therefore of information, 
apparition of shapes and quantization of these shapes (the cells can be considered as the 
'quanta' of life), spontaneous duplication and branching properties (see herebelow), etc... 
Moreover, due to the existence of a vacuum energy in quantum mechanics (i.e., of a non 
vanishing minimum energy for a given system), we expect the primordial structures to 
appear at a given non-zero energy, without any intermediate step. 

Such a possibility is supported by the symplectic formal structure of thermodynamics 
[110], in which the state equations are analogous to Hamilton- Jacobi equations. One can 
therefore contemplate the possibility of a future 'quantization' of thermodynamics, and 
then of the chemistry of solutions, leading to a new form of macroscopic quantum (bio)- 
chemistry, which would hold both for the prebiotic medium at the origin of life and for 
today's intracellular medium. 

In such a framework, the fundamental equation would be the equation of molecu- 
lar fractal geodesies, which could be transformed into a Schrodinger equation for wave 
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functions tp. This equation describes an universal tendency to make structures in terms 
of a probabihty density P for chemical products (constructed from the distribution of 
geodesies), given by the squared modulus of the wave function i/j = a/P x e*^. Each 
of the molecules being subjected to this probability (which therefore plays the role of a 
potentiality), it is proportional to the concentration c for a large number of molecules, 
P (X c but it also constrains the motion of individual molecules when they are in small 
number (this is similar to a particle- by-particle Young slit experiment). 

Finally, the Schrodinger equation may in its turn be transformed into a continuity and 
Euler hydrodynamic-like system (for the velocity V = {v+ + v^)/2 and the probability P) 
with a quantum potential depending on the concentration when P oc c, 

Q = -2I?2^. (106) 

This hydrodynamics-like system also implicitly contains as a sub-part a standard diffusion 
Fokker-Planck equation with diffusion coefficient V for the velocity It is therefore 
possible to generalize the standard classical approach of biochemistry which often makes 
use of fluid equations, with or without diffusion terms (see, e.g., [64, 119]). 

Under the point of view of this third representation, the spontaneous transformation 
of a classical system into a quantum-like system through the action of fractality and small 
time scale irreversibility manifests itself by the appearance of a quantum-type potential 
energy in addition to the standard classical energy balance. We therefore predict that 
biological systems must show an additional energy (quite similar to the missing energy of 
cosmology usually attributed to a never found 'dark matter') given by the above relation 
(106) in terms of concentrations, when their total measured energy balance is compared 
to the classically expected one. 

But we have also shown that the opposite of a quantum potential is a diffusion po- 
tential. Therefore, in case of simple reversal of the sign of this potential energy, the 
self-organization properties of this quantum-like behavior would be immediately turned, 
not only into a weakly organized classical system, but even into an increasing entropy 
diffusing and desorganized system. We tentatively suggest [96] that such a view may pro- 
vide a renewed way of approach to the understanding of tumors, which are characterized, 
among many other features, by both energy affinity and morphological desorganization. 

3.8 Duplication 

Secondly, the passage from the fundamental level to the first excited level now provides one 
with a (rough) model of duplication (see Figs. 13 and 14). Once again, the quantization 
implies that, in case of energy increase, the system will not increase its size, but will 
instead be lead to jump from a single structure to a binary structure, with no stable 
intermediate step between the two stationary solutions n = and n = 1. Moreover, if 
one comes back to the level of description of individual trajectories, one finds that from 
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Figure 13: Model of duplication. The stationary solutions of the Schrodinger equation in 
a 3D harmonic oscillator potential can take only discretized morphologies in correspondence 
with the quantized value of the energy. Provided the energy increases from the one-structure 
case (£"0 = 3T>U!), no stable solution can exist before it reaches the second quantized level at 
El = 5T>uj. The solutions of the time-dependent equation show that the system jumps from the 
one structure to the two-structure morphology. 



each point of the initial one body-structure there exist trajectories that go to the two final 
structures. In this framework, duplication is expected to be linked to a discretized and 
precisely fixed jump in energy. 

It is clear that, at this stage, such a model is extremely far from describing the com- 
plexity of a true cellular division, which it did not intend to do. Its interest is to be a 
generic and general model for a spontaneous duplication process of quantized structures, 
linked to energy jumps. Indeed, the jump from one to two probability peaks when going 
from the fundamental level to the first excited level is found in many different situations of 
which the harmonic oscillator case is only an example. Moreover, this duplication property 
is expected to be conserved under more elaborated versions of the description provided 
the asymptotic small scale behavior remains of constant fractal dimension Dp k, 2, such 
as, e.g., in cell wall-like models based on a locally increasing effective fractal dimension. 

3.9 Bifurcation, branching process 

Such a model can also be applied to a first rough description of a branching process 
(Fig. 14), e.g., in the case of a tree growth when the previous structure remains instead 
of disappearing as in cell duplication. 

Note finally that, although such a model is still clearly too rough to claim that it 
describes biological systems, it may already be improved by combining with it various 
other functional and morphological elements which have been obtained. Namely, one 
may apply the duplication or branching process to a system whose underlying scale laws 
(which condition the derivation of the generalized Schrodinger equation) include (i) the 
model of membrane through a fractal dimension that becomes variable with the distance 
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to a center; (ii) the model of multiple hierarchical levels of organization depending on 
'complexergy' (see herebelow). 




Figure 14: Model of branching and bifurcation. Successive solutions of the time-dependent 2D 
Schrodinger equation in an harmonic oscillator potential are plotted as isodensities. The energy 
varies from the fundamental level (n = 0) to the first excited level (n = 1), and as a consequence 
the system jumps from a one-structure to a two-structure morphology. 



3.10 Nature of first evolutionary leaps 

We have also suggested applications to biology of the new quantum-like mechanics in 
scale space [89]. 

In the fractal model of the tree of life described hereabove [19], we have voluntarily 
limited ourselves to an analysis of only the chronology of events (see Fig. 7), indepen- 
dently of the nature of the major evolutionary leaps. The suggestion of a quantum-type 
mechanics in scale space and of the new concept of complexergy [89, 101], which is a 
new conservative quantity appearing from the symmetry of the new scale variables (more 
precisely, of the fractal dimension become variable and considered as a fifth dimension) 
allows one to reconsider the question. 

One may indeed suggest that life evolution proceeds in terms of increasing quantized 
complexergy. This would account for the existence of punctuated evolution [41], and for 
the log-periodic behavior of the leap dates, which can be interpreted in terms of probability 
density of the events, P = l^/'p oc sin^[u;ln(T — Tc)]. Moreover, one may contemplate the 
possibility of an understanding of the nature of the events, even though in a rough way 
as a first step. 

Indeed, one can expect the first formation of a structure at the fundamental level 
(lowest complexergy), which is generally characterized by only one length-scale (this is the 
analog in scale space of the left part of Fig. 13 which concerns position space). Moreover, 
the most probable value for this scale of formation is predicted to be the 'middle' of the 
scale-space, since the problem is similar to that of a quantum particle in a box, with the 
logarithms of the minimum scale and maximum scale Ajv/ playing the roles of the walls 
of the box, so that the fundamental level solution has a peak at a scale \/Xm x Am- 
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The universal boundary conditions are the Planck-length /p in the microscopic domain 
and the cosmic scale L = A~^/^ given by the cosmological constant A in the macroscopic 
domain (see Sec. 3.1.2). From the predicted and now observed value of the cosmological 
constant, one finds L//p = 5.3 x 10^°, so that the mid scale is at 2.3 x 10^° /p ~ 40 /im. 
A quite similar result is obtained from the scale boundaries of living systems (~0.5 
Angstroms - 30 m). This scale of 40 //m is indeed a typical scale of living cells. Moreover, 
the first 'prokaryot' cells appeared about three Gyrs ago had only one hierarchy level (no 
nucleus) . 

In this framework, a further increase of complexergy can occur only in a quantized 
way. The second level describes a system with two levels of organization, in agreement 
with the second step of evolution leading to eukaryots about 1.7 Gyrs ago (second event 
in Fig. 7). One expects (in this very simplified model), that the scale of nuclei be smaller 
than the scale of prokaryots, itself smaller than the scale of eucaryots: this is indeed what 
is observed. 

The following expected major evolutionary leap is a three organization level system, 
in agreement with the apparition of multicellular forms (animals, plants and fungi) about 
1 Gyr ago (third event in Fig. 7). It is also predicted that the multicellular stage can 
be built only from eukaryots, in agreement with what is observed. Namely, the cells of 
multicellulars do have nuclei; more generally, evolved organisms keep in their internal 
structure the organization levels of the preceeding stages. 

The following major leaps correspond to more complicated structures then more com- 
plex functions (supporting structures such as exoskeletons, tetrapody, homeothermy, 
viviparity), but they are still characterized by fundamental changes in the number of 
organization levels. Moreover, the first steps in the above model are based on spheri- 
cal symmetry, but this symmetry is naturaly broken at scales larger than 40 /im, since 
this is also the scale beyond which the gravitational force becomes larger than the van 
der Waals force. One therefore expects the evolutionary leaps that follow the appari- 
tion of multicellular systems to lead to more complicated structures, such as those of 
the Precambrian- Cambrian radiation, than can no longer be described by a single scale 
variable. 

3.11 Origin of the genetic code 

We therefore intend, in future works, to extend the model to more general symmetries, 
boundary conditions and constraints. We also emphasize once again that such an approach 
does not dismiss the role and the importance of the genetic code in biology. On the 
contrary, we hope that it may help understanding its origin and its evolution. 

Indeed, we have suggested that the various biological morphologies and functions 
are solutions of macroscopic Schrodinger-type equations, whose solutions are quantized 
according to integer numbers that represent the various conservative quantities of the 
system. Among these quantities, one expects to recover the basic physical ones, such as 
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energy, momentum, electric charge, etc... But one may also contemplate the possibility 
of the existence of prime integrals (conservative quantities) which would be specific of 
biology (or particularly relevant to biology), among which we have suggested the new 
concept of complexergy, but also new scale 'charges' finding their origin in the internal 
scale symmetries of the biological systems. 

The quantization of these various quantities means that any such system would be 
described by a set of integer numbers, so that one may tentatively suggest that only 
these numbers, instead of a full continuous and detailed information, would have to be 
included in the genetic code. In this case the process of genetic code reading, protein 
synthesis, etc... would be a kind of 'analogic solutioner' of Schrodinger equation, leading 
to the final morphologies and functions. Such a view also offers a new line of research 
toward understanding the apparition of the code, namely, the transformation of what was 
a purely chemical process into a support of information and of its implementation, thanks 
to the quantization of the exchanges of energy and other conservative quantities. 

We intend to develop this approach in future works, in particular by including the scale 
relativity tools and methods in a system biology framework allowing multiscale integration 
[6, 101], in agreement with Noble's 'biological relativity' [65] according to which there is 
no privileged scale in living systems. 

4 Conclusion 

The theory of scale relativity relies on the postulate that the fundamental laws that govern 
the various physical, biological and other phenomenons find their origin in first principles. 
In continuity with previous theories of relativity, it considers that the most fundamental 
of these principles is the principle of relativity itself. The extraordinary success due to 
the application of this principle, since now four centuries, to position, orientation, motion 
(and therefore to gravitation) is well known. 

But, during the last decades, the various sciences have been faced to an ever increasing 
number of new unsolved problems, of which many are linked to questions of scales. It 
therefore seemed natural, in order to deal with these problems at a fundamental and first 
principle level, to extend theories of relativity by including the scale in the very definition 
of the coordinate system, then to account for these scale transformations in a relativistic 
way. 

We have attempted to give in this article a summarized discussion of the various de- 
velopments of the theory and of its applications. The aim of this theory is to describe 
space-time as a continuous manifold without making the hypothesis of differentiability, 
and to physically constrain its possible geometry by the principle of relativity, both of 
motion and of scale. This is effectively made by using the physical principles that directly 
derive from it, namely, the covariance, equivalence and geodesic principles. These prin- 
ciples lead in their turn to the construction of covariant derivatives, and finally to the 
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writing, in terms of these covariant derivatives, of the motion equations under the form 
of free-hke geodesic equations. Such an attempt is therefore a natural extension of gen- 
eral relativity, since the two-times differentiable continuous manifolds of Einstein's theory, 
that are constrained by the principle of relativity of motion, are particular sub-cases of 
the new geometry in construction. 

Now, giving up the differentiability hypothesis involves an extremely large number 
of new possible structures to be investigated and described. In view of the immensity 
of the task, we have chosen to proceed by steps, using presently known physics as a 
guide. Such an approach is rendered possible by the result according to which the small 
scale structure which manifest the nondifferentiability are smoothed out beyond some 
relative transitions toward the large scales. One therefore recovers the standard classical 
differentiable theory as a large scale approximation of this generalized approach. But one 
also obtains a new geometric theory which allows one to understand quantum mechanics 
as a manifestation of an underlying nondifferentiable and fractal geometry, and finally to 
suggest generalizations of it and new domains of application for these generalizations. 

Now the difficulty with theories of relativity is that they are meta-theories rather than 
theories of some particular systems. Hence, after the construction of special relativity 
of motion at the beginning of the twentieth century, the whole of physics needed to be 
rendered 'relativistic' (from the viewpoint of motion), a task that is not yet fully achieved. 
The same is true as regards the program of constructing a fully scale-relativistic science. 
Whatever be the already obtained successes, the task remains huge, in particular when 
one realizes that it is no longer only physics that is concerned, but now many other sci- 
ences, in particular biology. Its ability to go beyond the frontiers between sciences may be 
one of the main interests of the scale relativity theory, opening the hope of a refoundation 
on mathematical principles and on predictive differential equations of a 'philosophy of 
nature' in which physics would no longer be separated from other sciences. 
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